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I. INTRODUCTION 



For more than a decade various models on the behavior of quarks and gluons at small x has been confronted with 
a large amount of experimental data from HERA on the deep- inelastic scattering (DIS) structure function F2 [1-14]. 
In the small x regime, non-perturbative effects are expected to give a substantial contribution to F2. However, what 
is observed up to very low Q 2 ~ 1 GeV 2 values, traditionally explained by soft processes, is described reasonably well 
by perturbative QCD evolution (see for example [15]). Thus, it is important to identify the kinematical region where 
the well-established perturbative QCD formalism can be safely applied. 

At small x the Q 2 dependence of quarks and gluons is usually obtained from the numerical solution of the Dokshitzcr- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [16-20] [129]. The x profile of partons at some initial Qq and the 
QCD energy scale A are determined from a fit to experimental data [21-34]. 

On the other hand, when analyzing exclusively the small x region, a much simpler analysis can be done by using 
some of the existing analytical approaches of DGLAP equations in the small x limit [35-44]. In Refs. [35-37, 43, 44] it 
was pointed out that HERA small x data can be interpreted in terms of the so called doubled asymptotic scaling (DAS) 
phenomenon related to the asymptotic behavior of the DGLAP evolution discovered many years ago in [16, 17, 45]. 

In the present work we incorporate the contribution from higher twist (HT) operators of the Wilson operator product 
expansion to our previous analysis [44]. The semi-analytical solution of DGLAP equations obtained in Ref. [44] using 
a flat initial condition, is the next-to-leading order (NLO) extension of previous studies performed at the leading 
order (LO) in perturbative QCD [35, 43]. The flat initial conditions at some initial value Qq correspond to the case 
of parton distributions tending to some constant when x — > 0. 

In Ref. [44], both the gluon and quark singlet densities are presented in terms of the diagonal and ' components 
obtained from the DGLAP equations in the Mellin moment space. The '— ' components are constants at small x for 
any values of Q 2 , whereas the '+' components grow for Q 2 > Q\ as [130] 



exp 
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where a + = 4Ca/Po and b+ — 8[23CU — 26CV]Tr//(9/3o). In Eq. (1) and hereafter we use the notation a s = a s /(4n). 

The hrst two coefficients of the QCD /3-function in the MS-scheme are (3 = (11/3)CU - (4/3)Tr/ and /?i = 
(2/3)[17C^ — IOCaT^r/ — 6C^ Tr/] where / is the number of active flavors. This new presentation as a function of 
the SU(N) group Casimirs, with / active flavors, C A = N,T R = 1/2, T F = T R f and C F = (N 2 - 1)/(2N) permits 
one to apply our results to, for example, the popular N = 1 supersymmetric model. Of course, for TV = 3 one obtains 
the QCD result [44]. 

The analysis performed in our previous work [44] has shown very good agreement with HI and ZEUS 1994 data 
[4, 11, 12] at Q 2 > 1.5 GeV 2 . Here, we add the contribution from higher twist operators with the hope to describe 
also more modern data [2, 3, 7-10] at lower Q 2 

Moreover, in comparison with Ref. [44], in the present work we have solved the technical problem of "backward" 
evolution that leads us now to have the normalization scale Qq of DGLAP evolution in the middle point of the Q 2 
range. 



A. Basic formulae 



At this point of the introduction, we find convenient to present the basic results of our article: the twist-four and 
twist-six corrections to F 2 in the DAS approach. Thus, a reader who has interest only in application of the formulas 
to the analysis of F 2 can skip the following sections and start to read Section X, where the fits of F 2 are performed. 
We note, however, that some of the sections that follows contain also the contribution of power corrections to the 
derivatives dF 2 /d\nQ 2 and dlnF 2 /dln(l/x) and to the parton distributions (PD) (sec the Sections VI, VII, VIII 
and IX, respectively). 

The basic results of the present article are the twist-four and twist-six corrections to F 2 

F 2 (x,Q 2 ) = F^(x,Q 2 ) + ± F ;\z,Q 2 ) + -L F ; e (z,Q 2 ) , (2) 

where for the higher twist parts F% ' 6 BFKL-motivated evaluations [46-49] (in this case only the twist-four correction 
has been estimated) and the calculations [50] in the framework of the renormalon model (hereafter marked with 
superindex R) have been used. 

The latter case is essentially more complete and the predicted HT corrections can be expressed through the twist-two 
ones as follows 

F^\z,Q 2 ) = e ]T a?ft 4 (z,Q 2 )®/ T2 (*,Q 2 ) = E ( 3 ) 

a—q.G a—q,G 

where the symbol <g> marks the Mellin convolution (sec Eq. (55) below), the functions /x^ 4 (z,Q 2 ) are given in [50] 
and e = (X)i e f)// i s the average charge square for / active quarks. We call F 2 R ^ 4 and F^ 4 , the HT corrections 
proportional to the twist-two quark and gluon densities, respectively. 

Note that the parton distributions f^ 2 (z,Q 2 ) are multiplied by z, i.e., fg 2 (z,Q 2 ) — zq(z,Q 2 ) and f(?(z,Q 2 ) = 
zG(z,Q 2 ). Note also that we neglect the non-singlet quark density /a(z,Q 2 ) and the valent part fv(z,Q 2 ) of the 
singlet quark distributions, because they have the following small- a; asymptotics: /a (z, Q 2 ) <~ fv(z, Q 2 ) ~ x Xv , where 
Ay ~ 0.3 -j- 0.5. Thus, our quark density fa 2 (z, Q 2 ) contains only the sea part fs{z, Q 2 ), i.e. /„ 2 (z, Q 2 ) = fs{z, Q 2 )- 

For the leading twist part we have [44] at the LO and NLO approximations, respectively, 



Filoi^Q 2 ) = e.r q l {z,Q 2 ) , (4a) 
F 2 T2 (z,Q 2 ) = e (n 2 (z,Q 2 ) + ^La s (Q 2 ).r G 2 (z,Q 2 )) . (4b) 



Let us keep the NLO relation (4b) beyond the leading twist approximation. Then for the total F 2 (see Eq. (2)) we 
obtain 

F 2 (z, Q 2 ) = e (f q (z,Q 2 ) + ^la s (Q 2 )f G (z,Q 2 )] , (5) 



where f a (z, Q 2 ) are the parton distributions containing both the twist-two part [44] (see next Section) and the twist- 
four and twist-six contributions 
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For the HT part fa T4 ' 6 ( z > Q 2 ) calculations in the framework of the renormalon model have been used [131]. 

We would like to note that each HT term f* Tifi {z,Q 2 ) can be chosen in a quite arbitrary form and only the 
combination 

f q RTi ' e ^Q 2 ) + ^« s (Q 2 )/g t4 ' 6 (^0 2 ) (7) 

is unique, because we kept the original twist-two relation, Eq. (4b), to be same when HT corrections are incorporated 
(see Eq. (5)). 

Note that in our previous studies [51-53] we did not use the Eq. (5) to parameterize the HT corrections to F 2 . 
Instead we consider the following representation 

Ft Afi {z,Q 2 ) = e/f 4 < 6 (z,Q 2 ), (8) 

coming from the LO relation (4a) between F 2 and parton distributions. The choice (8) looks quite natural because 
HT corrections have been obtained in [50] at the LO approximation. However, this choice is only useful to fit F 2 data 
and it has no interest to study the parton distributions themselves: note that the HT corrections to the gluon density 
are absent in Eq. (8). Indeed, in the calculation of F 2 at NLO one has to take a gluon density as in the r.h.s. of the 
Eq. (4b). So, one should take the condition 

/g T4 ' 6 (^Q 2 ) = 0, (9) 

which is not so natural. Moreover, the choice (8) and (9) leads to a quite complicated form for the HT corrections to 
the quark density: there are two independent contributions ~ A T q 2 and ~ Aq (see Refs. [51-53] and formulas therein). 

In the work we also study x and Q 2 dependences of 0F 2 /d\nQ 2 and <91ni*2/<91n(l/x), that force to define the 
parton densities in a proper way. So, we take another quite natural choice 

C 4 ' 6 M 2 ) = ^V q ifi {z,Q 2 )®f T q \z,Q 2 ) = \F^{z,Q 2 ) , (10a) 

i.e., the HT quark (gluon) part of F 2 relates only to the corresponding quark (gluon) twist-two density. 

Note once again that the choice (10) corresponds exactly to the Eq. (5), i.e. to the extension of the standard 
twist-two relation (4b) between F 2 and parton densities at the NLO formulas with the purpose to include the HT 
contributions. 

Note also that for both of the above parton density choices the DGLAP equation will be violated by the HT 
corrections (see Section VI and discussions therein). 



B. Higher twist terms in the renormalon model 

As it has been already noted above it is useful to split the parton distributions in two parts 

f a (z,Q 2 ) = f+(z,Q 2 ) + f~(z,Q 2 ) , (11) 

where the both '+' and ' components contain twist-two and HT parts. 

The two component representation follows directly form the exact solution of DGLAP equation in the Mcllin 
moment space at the leading twist approximation (see [44]). 

The twist-two contribution is presented below in the Section II and the twist-four and twist-six parts can be 
expressed through the twist-two one as follows (here for simplicity we restrict our consideration by LO approximation): 

for the (singlet) quark distribution 

fq TA,+ ( z , Q 2 ) 64CfTr/ 4 . _ , . ( 



f? Ti '-{z,Q 2 ) _ 6AC F T R f 




(12b) 
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for the gluon distribution 

/* t4 ' + (z,Q 2 ) _ 8 a G 4 \ 2 U^p) n jQ 2 \\ , „, , (Uc] 

fo%^) -M^)\—oi^) +ln yw\)) + oip ^> (12c) 

glM!) = A ln ( + iz) (12d) 

/c 2 lo(^0 2 ) 5/9g «.(Q 2 ) ln W G \^\) + ° (Z) ' 

where a^ 4 are the magnitudes which should be extracted from the fits of the experimental data. The variables 
z a = zexp\p(u a )], where p(v a ) = + v a ) — ^(va)] and v a are the powers of the x — > 1 asymptotics of the parton 
distributions, i. e. f a ~ (1 — x)" a at x — > 1. From the quark counting rules we know that z/ q w 3 and ~ 4. Then, we 
get p{v q ) w 11/6 and p(^g) ~ 25/12, and there derivatives v\ v q) ~ —49/36 and y'iyo) w —205/144 (see Appendix B 
for further details). 

The functions /„ in Eqs. (12a, 12c) are related to the modified Bessel function I v and to the Bessel function J v by: 

T (n\ - \ if<T 2 = * 2 >0, 

^ (fT) - ti"J v (S), if o 2 = -S 2 < . (13) 

and the a and p values are given in the Section II by Eqs. (20) and (23) at the LO and by Eqs. (26) NLO, respectively. 

Note that the upper (down) line in the r.h.s. of Eq. (13) corresponds to the solution of the DGLAP equation for 
the "direct" ( "backward" ) evolution in the DAS approximation. 

The twist-six part can be easy obtained from the corresponding twist-four one as 



/^ 6 (z,Q 2 ) = -^x 



/^ 4 (z, Q 2 ) with a- 4 - a; 6 , ln ( - ln ' Q 



VK 6 



(14) 



The article is organized as follows. In Section II we shortly review basic formulae of the solution of DGLAP 
equation at small x values with the flat initial conditions, given in [44] . We show the possibility to add the backward 
evolution to the formulae. In Sections III and IV we present the set of formulae for the derivation dF^jd InQ 2 and 
for the effective slope <91ni*2/<91n(l/a;). Section V contains our suggestions about the contributions of twist-four and 
twist-six operators of the Wilson operator product expansion. In Sections VI-IX we consider the contributions of the 
HT operators to parton distributions and to derivatives of F2 in the framework of the infrared renormalon model. 
Section X contains the fits of experimental data for F 2 , predictions for its derivatives and some discussions of the 
obtained results. In the Appendix A we present Mcllin moments of renormalon contributions, calculated in [50], and 
obtain their contributions to the PD corresponding moments. In the Appendix B we illustrate the method [54, 55] of 
replacing at small x the convolution of two functions by simple products. The method is used in the present work for 
the correct incorporation of renormalon-type contributions of higher twists terms into our formulae. The conclusions 
contains summary of the results and suggestions about other applications of the approach. 



II. THE CONTRIBUTION OF TWIST-TWO OPERATORS 



As in [44] , we will work with the small x asymptotic form of parton distributions in the framework of the DGLAP 
evolution equations starting at some Qq with the flat function: 

fa(Qo) = A a (hereafter a = q, G) , (15) 

where A a are unknown parameters that have to be determined from data. 

We shortly compile below the main results found in [44] at the LO NLO approximations. 



A. Leading order 

The small x asymptotic results for PD, fl\ (a = q, G) and F^\ at LO of perturbation theory and at twist-two 
in the operator product expansion have been found in Ref. [44] : 

F^loi^Q 2 ) - e r q 2 LO (z,Q 2 ) , (16a) 

/; 2 loM 2 ) - /; 2 lSm 2 ) + /;, 2 loM 2 ) • 
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After Mellin inversion of the explicit moment solution to DGLAP equations, the '+' and '— ' PD components are 
given by: 

f G %o(z,Q 2 )= (fa +^ A A ?) 7 oKo) e-*+W>«> + O (p LO ) , (16c) 

f%£(*> Q 2 ) = ^ (^g + ^K) Plo Ii(oto) e- d +W™ + O (plq) , (16d) 

fcio^Q 2 ) = < e- d -W™ +0(z) , (16e) 

/; 2 lo (*> Q 2 ) = < e"^ 1 )- + , (16f) 



where 



are the regular parts of d+ and d_ anomalous dimensions, respectively, in the limit n — > 1. [132] 
We define the variable 

a s (Qo) 
a s (Q 2 ) 



At LO, in terms of the QCD scale Alo, it has the form: 

The argument ctlo m the LO is given by [133] 

a L o = 2yfdG G s LO ln(z) , (20) 

where 

dec - -4Ca//3o (21) 

is the singular part when n — ► 1 of rfcG = TggMA^A)) being 7g^(^) the LO coefficient of the gluon-gluon anomalous 
dimension. 

The prescription for the backward evolution given by Eq. (13) is the result, in the more general case, of the following 
representation of the series which appear in the inverse Mellin transformation of the exact solution for PD moments, 
(see for example Eq. (6) in Ref. [44]), 

V - - r^ 2 I (2Vt) = \tr v/2 i T ^ 2 ^ if * ^ ' (221 

^ o k\r(k + v + 1) ' IA2Vt} ~ 141 i J,(2^), if t < o . (22) 

And finally, in Eq. (16d) 



Let us note, that 



/ d GG SLO ctlo ,„„s 



p- v I v {p) - ^ ln y (l/z) at Q 2 -> Q 2 . (24) 
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Next-to-leading order 



The small x behavior of the twist- 2 parton densities fl 2 (a = q, G) and of FJ 2 at the NLO approximation has been 
presented in our previous paper [44]. Here we give the result that can also be used for Q 2 below the initial condition 
point Qq (where partons have the flat form in x as Eq. (15)) 



F^(z,Q 2 ) = e \fI\z,Q 2 ) + 



4T R f 



a s {Q 2 ) r G \z,Q 



The 



and 



r a 2 (z,Q 2 ) = r a 2 ' + (z,Q 2 ) + fz 2 '-(z,Q 2 ). 

PD components in the equations above are: 



f G 2 - + {z,Q 2 ) = A+(Q 2 ,Q 2 )I (a) exp(-J+(l) S - D+(l)p) 



f?' + (z,Q 2 ) =A+(Q 2 ,Q 2 ) 



O(p) , 
20C A 



(1 - d + _(l)a s (Q 2 )) ph (a) + a s (Q 2 )I (a) 



xexp(-d+(l)«-5+(l)p) + O(p) , 
f: 2 '~(z,Q 2 ) = A-(Q 2 ,Q 2 )cM-d-(l)s-D_(l)p) + O(z), 

where D±(n) = d±±(n) - (/3i//3o)d±(n); p = a s (Q%) - a s (Q 2 ) and 



a = 2\j (d+s + D+p) ]n(z 



(d+s + D+p) 



ln(z) 



21n(l/z) 



(25a) 
(25b) 

(25c) 



(25d) 
(25e) 



(26) 



A + G {Q 2 > Ql) = [1 - df_(l)a s (Q 2 )] A G 2 + §j [1 " d G + (l)a s (Q 2 ) - d G _(l)a s (Q 2 )} A? 
Aq(Q 2 , Qq) = A t q — Aq(Q 2 , Qq) , 



A+(Q 2 ,Q 2 ) = 



2T R f 



3C A 



AT? + 



C F 



|t2 



C a q 
a 8 (Q 2 )A+(Q 2 ,Q 2 ) . 



(27a) 
(27b) 
(27c) 

(27d) 



The different singular and regular parts of anomalous dimensions appearing in Eqs. (25)-(26) have the form [134]: 

(28a) 



d ++ = (23C A 26C F ) , d\_ = 



9/3o 



d%_ = 



d++(l) 



3/3 



C 2 



(36C(3) + 33C(2) 



4C F ((2) 



86 



C A 



.C F 



1643 
12 

547 C 2 

~m CF + 3 CA~ 



26 



T R f- 



26C F 



1 - 2 



C F 



d__(l) 



+ 



AC A C F 

Wo 
8C F 



C a 



1 - 2 



d* _(1) = C A 9 - 3-f - 4C(2) - - 1 - 2-f T R f , d%_(l) = 



.C f 



C A 
40C F T R f 



.Cf 
C a 



C a 

2C(3)-3C(2) + H+ 52I 1Z 2 



9C A 



3/3o 

dl + (l) = 0, d G + (l) 



T R f , 



C A 



1 - 2 



Cf 
C a 



T R f 



(28b) 
(28c) 

(28d) 
(28e) 



The corresponding numerical values are collected in Table I (see Ref. [44] for details) . 

We would like to note that the exact value of the variable o and the small x asymptotics of the modified Bessel 
function 



Iv(v) ~ exp (<r) at g — > oo 
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TABLE I: The values of the parameters used in the calculation of the parton distributions as a function of the number of flavors. 



/ 


d+ 


D+ 


d+(l) 


D+(l) 


d_(l) 


D_(l) 






d? + (l) 


3 


-4/3 


1180/81 


101/81 


-43.370269 


16/81 


1.974431 


2.779310 


80/27 


-29/9 


4 


-36/25 


91096/5625 


61/45 


-45.485532 


64/225 


3.108220 


2.618816 


320/81 


-89/27 


5 


-36/23 


84964/4761 


307/207 


-47.729779 


80/207 


4.674958 


2.458322 


400/81 


-91/27 


6 


-12/7 


8576/441 


103/63 


-50.057345 


32/63 


6.864360 


2.297828 


160/27 


-31/9 



are given in Introduction (see Eq. (1)) with \d + \ = a + and D + = b + + a+(/3i//3n). So, the most important part from 
the NLO corrections (i.e. the singlet part at x — > 0) is taken in a proper way: it comes directly into the argument of 
the Bessel functions and does not spoil the applicability of perturbation theory at low x values. 

We stress that the LO and NLO results given above coincide with the ones in Ref. [44] for positive values of s and 
slo (i- e - for the case Q 2 > Qq). 

Let us remind that these analytical expressions which have been obtained from the exact solution to the moment 
space DGLAP evolution equations in the asymptotic limit n — ► 1 have been already used in Ref. [44] to reproduce 
the small x behavior of parton distributions and lastly of DIS structure functions themselves. The consideration of 
negative values for s and slo leads us to apply the backward evolution in the present analysis and, thus, to have 
the possibility to choose any normalization point Q 2 and not only the low end of the Q 2 -evolution as it was done in 
Ref. [44]. 



III. THE CONTRIBUTION OF TWIST-TWO OPERATORS TO THE DERIVATIVE dF 2 /d\nQ 2 

In QCD the scaling violation of F 2 (z, Q 2 ) are caused by gluon brcmsstahlung from quarks and quark pair creation 
from gluons. In the low x domain the latter process dominates the scaling violations. F 2 is then largely determined 
by the sea quarks, whereas the dF 2 /d In Q 2 is dominated by the convolution of the splitting function P q Q and the 
gluon density. At the leading twist approximation the derivative dF 2 /d\nQ 2 relates strongly to the gluon distribution 
f G 2 (z, Q 2 ). Moreover, the derivative is measured with a good accuracy. Then, the dF 2 /d\nQ 2 experimental data can 
be successfully used to determine the characteristic properties of gluon distribution. 

The 8F 2 /d\nQ 2 data becomes even more important, when we add higher twist corrections into consideration. In 
the case of the twist-four terms (of sum of the twist-four and twist-six terms) in the renormalon model there are the 
two (four) additional parameters (see below Section IV) which may lead to problems to fit all them together only 
with help of F 2 experimental data. 



A. Leading order 



Note that at the LO approximation there are the following properties 



d 



dlnQ 2 
d 



—j— -ffc((TLo) 

Pto 



= 4CUa s (Q 2 ) -^rr /fc+i((TLo) , 



<91nQ 2 

which lead to the following results 

dlnQ 2 = asiQ) 
51nQ 2 - a ^> 



p£o/fcKo) = 4C A a s (Q 2 )p£o%-i|(<7Lo) (fc = 0, 1,2, . . .) , 



PLO lo(o-Lo) 



f G l + (^Q 2 )-Pod + (l)f£+(z,Q 2 ) 



df G \o^Q 2 ) (n2s $C F T R f , T 

dlnQ 2 = ^ 



3 

l T "' f fcioM 2 ) + O(z) 



+ 0( PL0 ) , 



iC A 



fq,Lo( Z iQ 2 ) (n2 , 8T R f , T 2,-/ n2 x , n( x 

— dln Q2 — = Ss W ) - — /g.lo( z ' Q ) + O(z) . 



(29a) 
(29b) 

(30a) 
(30b) 
(30c) 
(30d) 
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Thus, we have 

dF;l (z,Qi) 



dlnQ 2 



dlnQ 2 



ea s (Q 2 ) 



ST R f 



f£u>(z, Q z ) - 0o d+(l)/; 2 L +(z, Q 2 ) 



q,LO 1 



(31) 



The LO Q 2 evolution of the derivative dF^ 2 /dlnQ 2 is defined mostly by the corresponding evolution of the gluon 
distribution fa\o(z,Q 2 ), i.e. by the Eqs. (16b, 16c) and (16e). 



B. Next-to-leading order 



At the NLO approximation of perturbation theory the Eqs. (29) are replaced by 



d 



dlnQ 2 
d 



-M<t) 



dlnQ 2 

which leads to the following results 
= a s (Q ) 



a s (Q 2 ) [aC a ~ a s (Q 2 )f3od ++ \ -l_4 +1 (a) , 
a s (Q 2 ) Uc A -a s (Q 2 )/3 d ++ ] p k -%^\{<j) (k = 0, 1, 2, . . .) , 



(32a) 
(32b) 



dlnQ 2 



3C A 
80 C 



a: 



C a q 



AC A I {a) - p d+(l)ph(a) 



d++ + yC A (l + d~+(l)) 



U 2 '-{z,Q 2 ) 

dlnQ 2 



+ as{Q ^{f p h{a) ~ 
+ A, (4_(l)(l + «i + (l))-I ++ (l)) ph(a)} 



Al 2 (d_(l) + a s (Q 2 )d-(l)) 



- 4C A d g + _{l))l {a) 



exp(-J + (l)s - D+{l)p) + O(p) 



(33) 



= -(3 a s (Q 2 ) 



it2 



exp(-d_(l)s-D_(l)p) + O(z) . 



(34) 



Taking together equations (25), (30a), (30c), (33) and (34), after some algebra we have got the final result 



dFp(z,Q 2 ) 
dlnQ 2 



ea s (Q 2 ) 



8T R f 



(fh 2 (*,Q 2 ) + H^Q 2 )) - p d+(i)f T q 2 ' + (z,Q 2 ) 



a s (Q 2 )f3 d^(l)f^-(z,Q 2 ) 



where 



<P(z,Q 2 ) = <P+(z,Q 2 ) + $-(z,Q 2 ) , 
$+(z,0 2 ) = 0+(z,Q 2 ) exp(-d + (l) S - D + (l)p) + <D(p), 
$-{z,Q 2 ) = <f>-{z,Q 2 ) exp(-d_(l)*-£>_(l)p) + . 
The '+' and '— ' components in the equations above are: 

^{z,Q 2 ) =a s (Q 2 ) (Ag + ^A?} {f y W [d ++ + 2 -C A (l3 + 3d + (l)) 



(35) 



(36a) 
(36b) 
(36c) 



+ 4-d) 



- d°_(l))I (a) - 



A) 
±C A 



d\_{l) -d" ++ (l) ) p/iH }> + a s (Q 2 )A; 2 df + (l)/ (a) , (36d) 



(z,Q 2 ) = (a s (Q 2 ) a s (Q 2 )) |d?_(l) + - + 



^«,(0 2 )AI 2 



(36e) 



The values of the coefficients are given in Eqs. (28). 

Thus, the NLO Q 2 evolution of the derivative dF^ 2 /dlnQ 2 is defined mostly by the corresponding evolution of the 
gluon distribution fQ 2 (z,Q 2 ), i.e. by the Eqs. (25b, 25c) and (25e). 
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IV. THE CONTRIBUTION OF TWIST-TWO OPERATORS TO THE SLOPES OF F 2 AND OF 

PARTON DISTRIBUTIONS 

The behavior of F 2 and parton distributions can mimic a power law shape over a limited region of z,Q 2 : 



f a (z, Q 2 ) ~ z-^<5 2 ) and F 2 {z, Q 2 ) ~ z~ x ™^ . 

The slopes are effective ones because the parton distributions and F2 have mostly the Bessel-like form. 
Note that there are the following properties 



d 



d\n(l/z) 
d 

d\n{l/z) 



1 



T /fe_i(a) , 
p k+1 I k+1 (<j) (k = 0,1,2,...) , 



(37) 

(38a) 
(38b) 



which we will use below. 



A. Leading order 



The effective slopes have the form at the LO approximation 



•>eff,r2/ ^-,2^, 
A G,hO\ z ^ I 



fG%o( Z 'Q Z ) -^o(o'lo) 



— TT/'i" — 



A F2,hoi Z i Q ) 



PhO 



A(ctlo) 



(39a) 
(39b) 



The effective slopes A® ff and A*^, depend on the magnitudes A r2 of the initial PD and also on the chosen input 
values of Q\ and A. At quite large values of Q 2 >> Qq, where the '— ' component is not relevant, the dependence on 
the magnitudes of the initial PD disappear, having in this case for the asymptotic values: 



\ofi>2 / .~2\ ^l(cLo) 
A n t n *AZ, y ) = Plo ~ ~ PhO 



\eff,r2 / _ \cff.r2 / ri 2\ 



41n(l/z) ' 



^(ctlo) 
Plo ~- ~ Plo 

J 1 (^LO ) 



4 In (1/z) ' 



(40a) 
(40b) 



where symbol w marks approximations obtained by expansions of modified Bessel functions I n ((r). These approxi- 
mations should be correct only at very large o values (i.e. at very large Q 2 and/or very small x). It is the case (see 
Figs. 2 and 6). 

We would like to note that the slope Xp^LO a.s( z > Q 2 ) = A q S hOa.s( z i Q 2 ) coincides at very large a with one obtained in 
[57] (see also [15]) in the case of flat input. Note that the slope A^J^o as( z ' Q 2 ) 1& l &r g e then the slope A^' T L 2 as (z, Q 2 ) = 



x eff,r2 , n 2s x eff,T2 , n 2-, _ (h(<Jho) hi^ho) 

A G,hO,a,s\ Z > °» I ~ A F2,hO,as\ Z i ) ~ PhO 



Jo(olo) h(vho) 2 In (1/z) 



(41) 



that coincides with results of fits in Refs. [30, 34]. 



Next-to-leading order 



At the NLO approximation of perturbation theory we have the following properties of the effective slopes: the quark 
and gluon ones A® ff,r2 (z, Q 2 ) = d\nf£ 2 (z, Q 2 )/d\n (1/z) are reduced by the NLO terms that leads to the decreasing 
of the gluon distribution at small x. For the quark case it is not the case, because the normalization factor A T q 2,+ of 
the '+' component produces an additional contribution undamped as <~ (lnz) -1 . 
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Indeed, the effective slopes have the form, 

x ° (z ' 0) " r G 2 (^Q 2 ) P U*) ' (42a) 

A cff, T2(z Q2) _ / g r2 ' + (^Q 2 ) W (1 - 4-dK(Q 2 )) + (20C A /3) fls (Q 2 )7 1 (a)/p 2b 
q f?(z,Q 2 ) P I 1 (a){l-dl_(l)a s (Q 2 ))+(2QC A /3)a s (Q 2 )I (a)/p ' 

xeff> r2, ^ _ Af (z, Q 2 ) Q 2 ) + (4T fl //3) « S (Q 2 ) Xf(z, Q 2 ) f G 2 (z, Q 2 ) 

[Z ' Q) ~ f^ z ,Q 2 ) + (AT R f/Z)a s (Q 2 )f G 2 (z,Q 2 ) ' ^ 

The gluon effective slope A^ ' t2 (z, Q 2 ) is larger than the quark slope AJj ff ' r2 (z, Q 2 ), which is in excellent agreement 
with a recent MRS and GRV analysis [30, 34]. 
For the asymptotic values we have got 

\off,r2/ rs2\ ^1 (°0 1 r ao \ 

A G ,k s M) -Pj^ « P-^Jzy (43a) 

A e ff , T2(z Q 2 ) = W (1 - 4-(lMQ 2 )) + (20C A /3) Qs (Q 2 )7 1 (a)/p 
9 " S P / 1 ( ( r)(l-4_(l) as (Q 2 )) + (20C A /3)a s (Q 2 )7o( ( T)/p 

h(a) , 20C A 7 (q)7 2 (q) \ _ _J3_ lOCU a s (Q 2 ) 

7 2 (a) J ~ P 41n(l/z) + 3 pln(l/z) ' ^ j 

41n(l/z) ' 3 pln(l/z) " A ^ as 1 ' W > + pln(l/z) ' 1 j 

We would like to note that at the NLO approximation the slope A^'^ 2 (z, Q 2 ) lies between quark and gluon ones 
but closely to quark slope Af J 2 (z, Q 2 ), that is in agreement with Refs. [30, 34]. 
Indeed, 

>cff,r2, n 2N \cff>2, n 2, ( JlW , 26C A 2 \ / 7 (cr)7 2 (g) ^ 

V ' ' = ' ^) ~ 3~ J I f~~(a) — J 




1 , 26C A /o2 A 1 



G.as V-^W ) - ^F2,asV 

M l p - iinTTM + "T 1 ^ J 2pMiM ■ (44a) 

Q 2 ) - Af - 2 2 ) = 2C A a s( g 2 ) (l - - ■ (44b) 

Both slopes A^ ff ' r2 (z, Q 2 ) decrease with decreasing 0. A z dependence of the slope should not appear for a PD 
within a Regge type asymptotic (x _A ) and precise measurement of the slope A° ff,T2 (z, Q 2 ) may lead to the possibility 
to verify the type of small x asymptotics of parton distributions. The present data, however, are not enough to 
distinguish this slow x-dependence of A® ff,r2 (z, Q 2 ) (see Fig. 2). 

In the following Sections we study the higher- twist contributions to F 2 (x,Q 2 ), its derivatives and parton distribu- 
tions. 



V. THE HIGHER TWIST CONTRIBUTIONS FOR F 2 



In the Section we consider two different representations for twist-four effects. The first one comes from Regge-likc 
analysis [46-49]. Thus, it should have right asymptotics at x — > limit, but, unfortunately, the knowledge of its form 
is very restricted. 

The second one is based on the IR-renormalon model. The predictions can not reproduce the exact form of x — > 
asymptotics, calculated in Ref. [46-49] but gives rather good agreement with modern experimental data from HERA 
(see Section X). We think this agreement is similar to one (see Ref. [15]) at larger Q 2 values between DGLAP 
approach (even for its analytical simplification: the generalized DAS approach [44]) and experiment. 
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We would like to note here that in the analysis of experimental data performed below we consider both LO and 
NLO approximations in the twist-two case and for HT corrections in the rcnormalon case. In the BFKL-motivated 
approach, for simplicity [135] we restrict the calculation of the HT contribution to the consideration of LO Q 2 evolution 
alone. 

A. BFKL-motivated estimations for twist-four operators 

Twist-four operators are known [58] to have their own evolution equations but the diagonalization of the operator 
anomalous dimensions matrix is a very complicate problem. For our purpose, however, as the relevant limit is n — > 1, 
one can apply the results of Refs. [46-49], which have very simple form and are given in the classical DAS asymptotics 
considered in Section 2 of Ref. [44] 

Here we show that the contribution from twist-four operators can be represented in the same form as the twist-two 
operators by using the twist-four anomalous dimensions instead of the twist-two ones. 

For the singular part of twist-four anomalous dimensions we consider from Ref. [46] the result: 

lG %{n-l) = 2 7 ^((n-l)/2) (1 + e) , (45) 

where e is very small: e = 1/1224. 

Eq. (45) allows us to find the relation between the singular part of twist-four operators anomalous dimensions, 

7ob( n ) an d 7± 4 ("): w ith the twist-two ones, 7^ (n) and j±\n). It leads to the following relations: 

d T + 4 = d T G % = a 2 d+ = a 2 d GG , <T 4 = a 2 (L = , (46) 
where a 2 = 4(1 + e) and d + = d GG is given by Eq. (21). 

— t4 

The prediction for the regular parts d + (n) and d_(n) can not be obtained from Eq. (45), but it should be essentially 
less important in the kinematical range studied below. Then, in the analysis presented below, we proceed by fixing 
this non-singular part by means of a relation similar to Eq. (46): 

d\\\) = bd+(i), <r 4 (i) - 6d_(i), (47) 

and further we examine different "natural" choices of b: b = 0, 1 and a 2 /2. 

Note that the non-singular (when n — > 1) parts rf^ 4 (l), cT 4 (l) and d+(l), d-(l) determine the behavior of parton 
distributions and DIS structure functions at non-small x values. Usually fits to experimental data at intermediate 
and large values of x are performed with the help of the following forms for the structure function F 2 : 

F 2 (x,Q 2 ) = FZ 2 (x,Q 2 ) + ^F^(x) or (48) 

F 2 (x,Q 2 ) = FZ 2 {x,Q 2 )(l + ^/ 2 T4 W) (49) 

with Q 2 -independent functions FJ 4 (x) or /J 4 (x). 

In fact Eq. (48) is closed to our choice 6 = 0, i.e. the twist-4 contribution docs not evolve logarithmically with 
Q 2 . Also Eq. (49) is analogous to the choice 6=1, i.e. twist-two and twist-four operators have the same logarithmic 
Q 2 -dependence at large and intermediate x values. Lastly, the choice b = a 2 /2 corresponds to the hypothese about 
applicability of Eq. (45), obtained in the classical DAS limit, to a more wide generalized DAS approximation considered 
here. 

By analogy with Section II we represent the twist-four contribution split in the '+' and '— ' parts: 



Fl\z,Q 2 ) - e/; 4 (z,Q 2 ) , (50a) 

/; 4 (z,Q 2 ) = f; 4 > + (z,Q 2 ) + f: 4 '~(z,Q 2 ) . (50b) 
The '+' and '— ' PD components are: 

f G 4 - + (z, Q 2 ) = (a£ + %{aa w ) e^ 1 ^ + O (p LO ) , (50c) 

/; 4 ' + M 2 ) = Hf ( A o+^<) ^Lo7i(a<7Lo)e-^«-° + 0( PW ) , (50d) 

f^-(z,Q 2 ) = -^Al 4 e- bd -^o + o(z) , (50e) 

f T q 4 '-{z,Q 2 ) - A T q A c - bd -W^o + 0{z) , (50f) 
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because the corresponding twist-four projectors (see [59]) have the following form [136]: 

b 



p t4,+ t4,- = + Ji t4,- i _ t4,+ 

99 GG 99 2 ' aa aa ' 

fc (51) 
ri,± _ ± t4,± _ ± 

£ 9 G - £ 9G ^2" ' £ Gg ~ £ G 9 ■ 



In Eqs. (50c-50f) the twist-four parameters A T a A (a = q, G) have to be determined from fits to experimental data. 
The full contribution (i.e. the sum of twist-two and twist-four parts) is given by: 



U-.Q-) = /;; 2 (: Q-) + J:\z.Q 2 ) and (52) 

1 



F 2 (z,Q 2 ) = FJ 2 (z,Q 2 ) + —F^(z,Q 2 ), (53) 

where the leading twist contributions fa 2 (z,Q 2 ) and F 2 2 (z,Q 2 ) are given at LO by Eqs. (16a, 16b) and at NLO by 
Eqs. (25a,25b). 



B. Renormalon model predictions for twist-four operators 



The full small x asymptotic results for parton densities and F 2 structure function in the framework of the infrared 
renormalon model, i.e. F R , at LO of perturbation theory in the twist-four part: 



1 



F 2 R (z,Q 2 ) = FZ 2 (z, Q 2 ) + — F R ^(z,Q 2 ) , 



(54) 



where F^iz, Q 2 ) is given by Eqs. (16a-16f) at the LO approximation and by Eqs. (25a-25e) at NLO one, respectively. 
The twist-four term F 2 Rt4 (z,Q 2 ) has the form (3), i.e. 

- e F^\ Zl Q 2 ) = £ a :*ji:\z,Q 2 )®r a 2 (z,Q 2 ), 

a—q,G 

where the symbol <S> marks the Mcllin convolution 

A(z) ® B(z) = f ^ A{y) B (-\ . 

J z y \ v / 

The corresponding Mellin transforms of /x£ 4,6 (z, Q 2 ) 

p^ 6 (n,Q 2 ) = f dzz n - l K ifi (^Q 2 ) 
Jo 



(55) 



(56) 



are presented in the Appendix A (see Eqs. (A2-A4) and (A7)). 

Looking the n-space representations for renormalon power-like corrections given in Appendix A and applying the 
technique to transform the Mellin convolutions to standard products at small x (see [54, 55] and Appendix B) we can 
represent the Eq. (3) in the form 



1 



(z,Q 2 ) 



64T fl / 
15/3 2 



l I Y-l 101 1 , 
&r<8 1 + — + - In 
G ' 120 2 



01 

c.t4\ 



2C F aI 



2 11 - ! 2291 1 , 
-5 1 - -TT-r + - In 



' Q 2 \ 



120 



3600 



a r4 



5- 1 - 



139 
120 



s:\zm 2 ) 



The operators S 1 and 5 2 are defined as follows (see Appendix B for details) 

[r a 2 -(z,Q 2 )] = j-r a 2 >-(z,Q 2 ) , s- 2 [r a 2 >-(z,Q 2 )] = ^/; 2 >-(*,q 2 ) 



s- 1 



p k 1 I\ k - 1 \{a) 



p k I k (a) = P k - 2 l\ k - 2 \(o) . 



(57) 

(58a) 
(58b) 



13 



Note that the Eqs. (16) and Eqs. (25) have been obtained in [44] with the accuracy O(p) for the '+' component 
and with one 0{z) for the '— ' component, respectively. It leads to the fact that we should use only the most singular 
terms in the r.h.s. of Eq. (57): i.e. the terms (5 _1 and <~ ln((3 2 /|a G 4 |) for the gluon part and the terms 5~ 2 and 

ln(Q 2 /|a^ 4 |) <5 _1 for the quark part. 

Then, the Eq. (57) should be replaced by 



-f Rt Hz o 2 ) - MTlif 

e 2 [Z ' Q) - 15/3 2 



a r4 

a G 



+ 2C F a; 



r4 



-In 



(59) 



Applying the operators S 1 and S 2 separately to the '+' and '— ' components of fa 2 (z, Q 2 ), we obtain the following 
results for F 2 Rt4 (z, Q 2 ): 



where 



-(z,Q 2 ) + F 2 RTi '~(z,Q 2 ) , 



(60a) 



? Rt4,+ 



(z,Q 2 ) = 



+ 



32T R f 
15/3 2 



a G 4 



2 h{a) 
pTo(<?) 



In 



Q 2 



>t4| 



4W <(1-4_(1K(Q 2 )) 



- Fr 



RtA, 



+ 



~(z,Q 2 ) = 



3C A 9 v 
-^—a s (Q ) 



2Ha) + J Q 2 



32T fl / jr2 



15/3 2 



2C A sZ 



In 



In 



a G *ln 

( Q 2 



t4\ 
1 I 

)2 



Phil) 



+ ln 



r4 



z G |a G 



r4 



z <? | a <j 



Lt4| 



(60b) 



(60c) 



C. Incorporation of twist-six contributions in the framework of the renormalon model 



We shortly demonstrate the twist-six contributions in the framework of the renormalon model. 
When we added the twist-six part, the full small x asymptotic results for PD and F 2 en structure function at NLO 
of perturbation theory: 



F R (x,Q 2 ) = FZ 2 (x,Q 2 ) + ±F Rt4 (z,Q 2 ) + -L F ™( Z ,Q 2 ), 



By analogy with twist- four case the twist-six term f Rr6 (z, Q 2 ) has the form: 

-F R ^{z,Q 2 ) = ]T ^K 6 {z,Q 2 )®fl 2 {z,Q 2 ) , 



(61) 



(62) 



a—q,G 



where Ha 6 {z, Q 2 ) are given in [50]. The corresponding Mellin transform of /^ 6 (n, Q 2 ) is presented in the Appendix A 
(see Eqs. (A2), (A5,A6) and (A7)). 

By analogy with the previous subsection applying the technique to transform the Mellin convolutions to the standard 
products at small x (see [54, 55] and Appendix B), we can represent the Eq. (62) in the form 



a G 



~ i 2663 1 , 

S^ 1 H h - In 

3360 2 



2CW 



t6 



1 



143 - J 870637 
3360 ~ 1411200 + 2 n 




fa(z,Q 2 ) 



s- 1 - 



3217\ 



3360 



(63) 
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Considering only the most singular terms in the r.h.s. of (63), i.e. the terms S 1 and ~ ln(Q 2 / i/|ajj 6 |) for the gluon 
part and the terms S~ 2 and ln(Q 2 /y^|aJ 6 |) <5 _1 for the quark part, we obtain immediately the following results: 



7 15/? 2 



a G 6 j^ + iln(^)j /5 2 (z,Q 2 ) 



(64) 



which is very close to the twist-four one, see Eq. (59): 
1 



-F R ^(z,Q 2 ) = --x 



f Rr \z,Q 2 ) with a^a^ 6 , In 



Q 2 



In 



2 



(65) 



Note that the representation (65) of the twist-six terms in the terms of the twist-four ones is universal and has 
quite compact form and, thus, it will be often used below. 

Because the forms of the twist-four and twist-six contributions are very similar, it is possible to present quite 
compact form for the full contribution of the higher- twist operators F RhT (z, Q 2 ) 



where 



and 



F R (z 7 Q 2 ) = F^(z,Q 2 ) + F RhT { Zl Q 2 ) , 
F 2 Rh ^(z,Q 2 ) - F Rh ^ + (z,Q 2 ) + F Rh ^-(z,Q 2 ) 

a£ m \2h{a) 



(66a) 
(66b) 



\ F Rh ^{z,Q 2 ) =^.^ + (z,Q 2 ) £ k 



In 



Q(m-2) 1 p7o(f7 ) V|a™l 



Q 2 



+ AC f Rf Jf. ( (i - 4-dK(Q 2 )) \ 2 - J M + J Q2 



+ 



3C A Q(™-2) 
20C A 



I P 2 hi?) 



In 



Plo(cr) 
Q 2 \ hW) ' 



If r ^-(z,q 2 ) = 3 -^fr G 2 -( Z ,Q 2 ) E *- 



In 



P^o(o-) J 

Q 2 



(66c) 



Q(m-2) | a rm|P n 



2CV 



ln — In 



l Q(™- 2 ) I v -., 

where fc 4 = 1, = —8/7 and p 4 = 1, p$ = 1/2 



Q 2 



(66d) 



VI. THE HIGHER TWIST CONTRIBUTIONS FOR THE DERIVATIVE dF 2 /d\nQ 2 

By analogy with the previous Section we consider firstly only the twist-four terms in the framework of the infrared 
renormalon model. The contribution of the twist-six terms will be incorporated shortly at the end of this Section. 



A. Renormalon model predictions for twist-four operators 



Note that there are the following properties 
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where we keep only most important terms (see discussions in the previous Section and Eq. (59)). 
In this approximation we easily obtain that 



dF R (z,Q 2 ) _ 9FJ 2 (z,Q 2 



91nQ 2 



<91nQ 2 



1_ f dF R ^(z,Q 2 ) 
Q 2 \ dlnQ 2 



-Fp\z,Q 2 ) 



and 



dF^{z,Q 2 ) 



d\nQ 2 ~ 3 

The value of F 2 Rri (z, Q 2 ) is given by Eqs. (60a)-(60c) and 



e^la s (Q 2 )^ 4 (z,Q 2 ) 



*% T \z,Q 2 ) 



16CU „ t2i+ 



2 J G 



5/3 2 



AC F T R f A 2 I 2 (a) 



2 h(<j) 

f Q 2 



In 



a r4 



In 



(68a) 



(68b) 



(68c) 



Thus, we see that the twist-four corrections to F 2 and dF 2 /d\nQ 2 have opposite signs, because dF Rri / dhiQ 2 
a s (Q 2 ) and the most important twist-four contribution is given by F Rt4 (z 7 Q 2 ). 



B. Incorporation of twist-six contributions in the framework of the renormalon model 



Following to the subsection V C of the previous Section and considering the properties 



1 



1 



•In 



Q 4 \Q 



A 



2 \ 1 



d\nQ 2 Q 4 Q 4 ' dlnQ 2 

together with the one (29), we immediately obtain that 
dF R (z,Q 2 ) dFp(z,Q 2 ) , 1 fdF™(z,Q 2 ) 



1 

Q 2 



2 In 



In 



Q 2 \Q 



A 5 



<91nQ 2 



and 



<91nQ 2 Q 2 V dlnQ 2 

dF R ^(z,Q 2 ) _ _ ^8T fl / 



1 fdF™(z,Q 2 ) 



Q 4 V 91nQ 2 



<91nQ 2 

The value of .f flr6 (z, Q 2 ) is given by Eq. (65) and 



= e ■ 



as (Q 2 )$^(z,Q 2 ) . 



$g rt (z,Q 2 ) = - -x 



$^ 4 (z,Q 2 ) with aj 4 ->a?>, In 



«1 

| a r4| 



hi 



^T 6 



(69) 



2F 2 ^ b (z,Q 2 ) (70a) 



(70b) 



(70c) 



Thus, we see that by analogy with the case of F 2 (z, Q 2 ) itself, for the derivation (70) the twist-six terms partially 
compensate the contributions of the twist-four terms. 



VII. PARTON DISTRIBUTION FUNCTIONS IN THE RENORMALON MODEL APPROACH 

It is clearly to see that the standard parton distributions f q (x,Q 2 ) and fa{ x , Q 2 ) fitted with help of experimental 
data do not coincide with the above twist-two ones fq 2 (z, Q 2 ) and Jq 2 (z, Q 2 ). These PD f q (z, Q 2 ) and fc(z, Q 2 ) arc 
usually defined keeping their twist-two relations (16a) or (25a) with the structure function F 2 (z,Q 2 ), i.e. 

At LO 

F 2 {z,Q 2 ) = ef q (z,Q 2 ) , (71) 
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At NLO 

F2(z,Q 2 ) = e (f q (z,Q 2 ) + ^la s (Q 2 )f G (z,Q 2 )) . (72) 



Thus, the parton distributions f q (z, Q 2 ) and /g(z, Q 2 ) can be strongly deviated for the corresponding the twist-two 
densities fg 2 (z, Q 2 ) and fQ 2 (z 7 Q 2 ) at quite low Q 2 values, because there are the HT corrections to F^iz, Q 2 ). 

The HT correction to the parton dstribution at the LO was presented in the Introduction already. Here we present 
the results at the NLO. As it was in the previous Section, we consider firstly the twist-four corrections. 

A. Twist-four corrections to (singlet) quark distribution 

Consider firstly the (singlet) quark parton distribution f g (z,Q 2 ). From the Eq. (16a) and the analysis of the 
Section V we can obtain that 

f?(z,Q 2 ) = f; 2 (z,Q 2 ) + ^f^\z,Q 2 ) , (73a) 



where f^ Ti {z,Q 2 ) is given at the LO by Eqs. (12a) and (12b). 

omplcte expressions direct] 

jf(z,Q 2 ) = if+(z,Q 2 ) + ff-frQ*), (781.) 



It is useful to represent also the complete expressions directly for f?(z, Q 2 ). 



where at the NLO 




+ 0{p) , (73c) 



f^ + {z,Q 2 ) = i | QAC F T R f a^ 4 j 2 h(a) (l - 4_(l)a s (Q 2 )) + (20C A /3)a s (Q 2 )I 2 (a)/ p 
f q 2 -' + {z,Q 2 ) ' + 15/3 2 Q 2 \p 2 T 1 (a)(l-d q + _(l)a s (Q 2 ))+(20C A /3)a s (Q 2 )I (a)/p 

1 IqQt) (1 - 4_(l)a s (Q 2 )) + (20C A /3)a s (Q 2 )I 1 (<j)/p } 
Ph(a) (1 - d\_{l)a s (Q 2 )) + (20C A /3)a s (Q 2 )I (cr)/p J 

f?'~M 2 ) , . 'iifyvw a-' J. /IV ( Q 2 , , ■ 

/^■-( Z Q2) =1+ lr ^ 2 7S7) ln l7) ln l-L7jl-PW[ + O(z) . (,,-!<!) 

We clearly see that the twist-four terms are responsible for the additional positive contributions to the quark 
distribution, which are very important at low Q 2 values. 

So, the experimentally extracted quark distribution f q (z, Q 2 ), which have the leading twist relations (71) and (72) 
with F 2 (z, Q 2 ), strongly deviates form the leading twist quark distribution fg 2 (z, Q 2 ). At quite low Q 2 values, where 
fq 2 (z,Q 2 ) had the quite flat behavior closed to (15), the full quark distribution f^(z,Q 2 ) will rise at z —* (see 
Eqs. (73c) and (73d)), because a^ 4 > (see Tables VI ,VII). This rise is in full agreement with the corresponding 
experimental data (see Tables VI ,VII, Figure 9, Section X and discussions therein). 

B. Twist-four corrections to gluon distribution 

Consider now the gluon parton distribution fc{z,Q 2 )- From the Eq. (16a) and the analysis of the Section V we 
can obtain that 

f§(^Q 2 ) = fai^Q 2 ) + ^f£ T \z,Q 2 ) ■ (74a) 

where J* r4 (z, Q 2 ) is given at the LO by Eqs. (12c) and (12d). 

For the gluon distribution in the NLO we have the similar relations 

fBi^Q 2 ) = /g' + M 2 ) + Q 2 ) , (74b) 

IM = 1 + 1 a " JgW +ln /QM| +0{) (74c) 

f G 2 ' + (z,Q 2 ) + 5/3 2 a s (Q 2 )Q 2 1 ' " ' f ' 



fa'~(^Q 2 ) W 2 a s (Q 2 )Q 2 \z%\^ 




^ (Z ' Q2) =l+iT77^>lnf-^) - Cm:) . ,74cl, 
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So, as in the case of the quark distribution, the experimentally extracted gluon density fa(z,Q 2 ), which has the 
leading twist relation with F 2 (z,Q 2 ) and dF 2 /d\nQ 2 , strongly deviates form the leading twist gluon distribution 
fc?{ z iQ 2 )- At quite low Q 2 values: Q 2 <~ Qq , where fa 2 (z,Q 2 ) had the quite flat behavior closed to (15), the full 
gluon distribution f^(z, Q 2 ) falls at x — > 0, because sJq < (see Tables VI ,VII). The behavior is in full agreement 
with the corresponding experimental data (see Tables VI ,VII, Figure 9, Section X and discussions therein). 



C. Twist-six corrections to parton distributions 



We shortly demonstrate the twist-six contributions to parton distribution in the framework of the renormalon 
model. When we added the twist-six part, the full small x asymptotic results for parton distributions is 

fa(z,Q 2 ) = f T a 2 (z,Q 2 ) + l^f^\z,Q 2 ) + if /« T6 (s, Q 2 ) = f T a 2 (z,Q 2 ) + f^(z,Q 2 ) , (75) 
where fa TQ {z,Q 2 ) are given by Eqs. (14): 



f a ^(z,Q 2 ) = --x 



fa \ z , Q ) with a^-af , In 



Q 2 



In 



Q 2 



The twist-six corrections do not change the results for parton distributions obtained in the previous subsection. 

VIII. THE HIGHER TWIST CONTRIBUTIONS TO THE SLOPES OF F 2 AND OF PARTON 

DISTRIBUTIONS 



Consider the power-like corrections to the twist-two effective slopes X C p 2 T2 (z,Q 2 ) and A° fi,r2 (z, Q 2 ) (a — q,G) 



introduced in the Section IV. The effective slopes have the following form 

d 



^ C F2{ z tQ 2 ) — 

v, tt (z,Q 2 ) = 



d\n{l/z) 

• n*\ = - 

d\n{l/z) 



F T 2 2 {z,Q 2 ) + 1- 2 F 2 R ^Q 2 ) + ^Ft\z,Q 2 ) 

f r a 2 (z,Q 2 ) + ^L RT \Z,Q 2 ) + ±f™{z,<?) 



(76) 
(77) 



Using Eqs. (38), the derivations dF 2 2 / d\n{\/ z), df^ 2 /<91n(l/ ' z) and 9/ ( f Tm /91n(l/z), (to = 4,6), can be repre- 
sented as the sum of two components ('+' and '— ') which are obtained from the corresponding ('+' and '— ') PD 
functions. One can show that 

df^+(z,Q 2 ) _ MC F T R f T4 J 2 Jo(tr) (l - d q + _(l)a s (Q 2 )) + (20C A /3)a s (Q 2 )h(a)/ p 



dln(l/z) 



15/3 2 q [Ph(a) (1 - rfl_(l)a s (g 2 )) + (20C A /3)a a (Q*)I (a)/p 



+ ln 



,t4 



df^~{z,Q 2 ) MC F T R f . 

a' In 



Q 2 



<91n(l/z) 

9/g T4 ' + (z,Q 2 ) 

01n(l/z) 

dft 4 -(z,Q 2 ) 

<91n(l/z) 

dfa T6± ( z ,Q 2 ) 

<91n(l/z) 



15/3 2 



q I q \ I 



5/3 2 «4Q 2 ) 



p!M\r G 2 '+( Z ,Q 2 ), 



a, 



t4\ 



jj a G 4 /t2,-/ n 2N 

5/? 2 a s (Q 2 ) /G 

dfa Ti,± (z,Q 2 ) 



-7 X 



d\n{l/z) 



with a 



r4 1 a^ 6 , In . 



Q 2 



In 



g 2 



(78a) 
(78b) 
(78c) 
(78d) 
(78e) 



The Eqs. (76) and (77) together with the Eqs. (42) and (78) give a complete information about the full and 
asymptotical values of the slopes \p\(z, Q 2 ) and A° ff (z, Q 2 ). The results will be demonstrated on Figs. 2, 6 and 7. 
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It is possible, however, to give a simple demonstration of the effect of the HT corrections. Following the Sec- 
tion IV, we can prepare also the results for the higher twist corrections to the asymptotical values of A^ 2 ' r2 (z,(5 2 ) 
and A° ff ' T2 (z, Q 2 ), which can be obtained by neglecting the '— ' components. Restricting ourselves by the twist-four 
case we can estimate the value of the slopes as (z, Q 2 ) and X^ as (z 7 Q 2 ) in the form 



^ C F2,hs( z tQ 2 ) — ^F2,as( Z ' Q 2 ) 



Q 



,c3,Rt4 
2 A F2,as 



(z,Q 2 ) + 



1 

Q 1 



where at the LO 



x eflF,flr4/ n 2\ 
A F2,as l z > J 



16CU 
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5/? 2 2p L0 ln(l/z) 




+ ln 



\cff,i?r4 



( z ,Q 2 ) = 



q 



.^o(CTLo) - ^2(o"LO 
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64C F T R f a; 
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A 



cff,i?r4 
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15/3g 2p LO Ml/z) \ p LO 
5/3 2 a,(Q 2 ) ^ I^lo), 
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(79) 
(80) 



AVlo) 
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+ ln ^7 



5/3 2 a s (g 2 ) /9 L0 ln(l/^) 



(81a) 
(81b) 
(81c) 
(81d) 
(81e) 



From the equations (81b) and (81d) it possible to see that the slopes A^ as (2;, Q 2 ) and Xq S as (z, Q 2 ) have got the 
positive twist-four corrections, that is in full agreement with the corresponding experimental HI and ZEUS data for 
the slope A_F2 at low Q 2 values (see Fig. 7). However, the difference between the twist-four corrections to these slopes 
is negative , because ag 4 < (see Tables 4, 6-8): 



>cff,i?T4/ n 2\ 
A F2,as \ Z i *°i I 



X 



off,i?r4 
q,as 



(z,Q 2 




5p 2 pLoln(l/2) 



(82) 



Thus, the inequality Xp\ as (z, Q 2 ) > Xq S as (z, Q 2 ) coming form Eq. (44b) takes place for not very small Q 2 , because 
it is suppressed by power corrections. 

We would like to note that the equations (81) are valid only at not very small Q 2 values, where we can neglect 
the terms ~ 1/Q coming from expanding the denominator and from the twist-six terms. The small Q 2 behavior of 
Ajf^-z, Q 2 ) can be easy demonstrated at the point Q 2 — Q 2 in the following Section. 



IX. PARTON DISTRIBUTIONS IN THE RENORMALON MODEL AT Ql 

As it has been already shown in the previous Section the total PD functions f q {z,Q 2 ) and fciz^Q 2 ) fitted in 
experiments data do not coincide with the above twist-two ones fq 2 (z,Q 2 ) and j^{z,Q 2 ). It is very useful to 
demonstrate the difference at Q§, at the starting point of the DGLAP evolution. 

We begin the analysis with the consideration only the twist-four terms. The results can be calculated from the final 
formulae of the previous Section but it is simpler to repeat all calculations given in the Section V. At Q 2 — Qq all 
results simplify essentially because the leading-twist parton distributions are constant A q and A G at the point. 
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A. Parton distributions at Qq 

From the Eqs. (15) and (73) we can easy obtain at Q 2 = Q 2 ,, that 



1 



1 



f a (z,Ql) = Al 2 + f^( z ,Ql) + /flre^gg) 

1^0 (^o 



where at the LO 
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(83a) 



(83b) 
(83c) 
(83d) 



Thus, the total parton distributions f q (z, Ql) and fa{ z , Qo) are strongly deviated for the corresponding the twist- 
two densities /^" 2 (z, Q§) = A q 2 and f G 2 (z, Ql) = ^g 2 - Because usually the fitted values of & q 4 (a G 4 ) are positive 
(negative), the twist-four terms lead to positive and negative contributions in the case of quark and gluon densities, 
respectively. The twist-six terms do not change the results essentially. 



B. The effective slopes of F2 and of parton distributions at Qq 

To estimate the values of the effective slopes at low Q 2 values we can look on their behavior at Qq, where our 
formulae simplifies essentially. In the approximation, when the twist six contributions arc negligible, we can easy 
obtain from the Eqs. (83) 



\^' R (Z, Qq) 

X S p2 R (z, Qq) 
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a G G 



C F Af 



(84a) 
(84b) 
(84c) 



Because a G 4 < 0, it is easy to see that A^'^z, Qq) < X q S ' R (z, Qq). This indicates that the inequality Xp\(z, Q 2 ) > 
X q (z,Q 2 ) seems to be correct only at quite large Q 2 values (see also the previous section and discussions therein), 
where the twist-two terms give basic contributions. 

Note also, that at Qq the slope X q S ' R (z, Qq) rises at x — ► 0, but the gluon slope Xq' R (z,Qq) is negative and x- 
independent. Thus, Xf' R (z,Q 2 ) > Xf' R (z,Ql) at low Q 2 that is in full agreement with the recent experimental data 
from HERA (see, for example, the review [15]). The twist-six terms do not change the above results essentially. 



X. RESULTS OF THE FITS 



With the help of the results obtained in the previous sections we have analyzed F2(x, Q 2 ) HERA data at small x 
from the HI [1-6] and ZEUS [7-14] collaborations as separately, as well as together. 

Without higher-twist corrections our solution of the DGLAP equations depends on five parameters, i.e. Qq, Xq, 
Aq, A t2 and A (or, equally well, on a s (Mz))- The incorporation of twist-four and twist-six corrections leads to two 
and four additional parameters, respectively. 
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In order to keep the analysis as simple as possible we have fixed A^jg to the values given in Eq. (85), which 
corresponds to a s (Mz) — 0.1166, obtained recently by ZEUS [7]. The analyzed data region was restricted to x < 0.01 
to stay in the kinematic region where our results are expected to be applicable. The \ 2 minimizations were done 
with MINUIT [60]. In the fits the errors are statistical and systematical added in quadrature. Finally, the number of 
active flavors was fixed to / = 3 and 4 for comparison. 

A. Leading twist approximation 

Tables II and III summarize the results of the fits to HI and ZEUS data using twist-two formulas at LO (16) and 
NLO (25) approximations. 

We can see in Tables II and III and in Fig. 1 that the qualities of the fits are very similar for the LO and NLO 
approximations. This suggest that perturbation theory works well in the small x regime. This is in accord with 
Refs. [61-63] (see also recent review [64]), where it was shown, that the argument of the strong coupling constant is 
effectively much larger as Q 2 in the small x domain. 

Hovewer the similarity of the results found at LO and NLO fits does not agree with our previous analysis [44] , where 
NLO corrections essentially improved the comparison between QCD and experiment. This disagreement relates mostly 
to the incorrect use of the same value of the QCD parameter A in [44] in both LO and NLO cases. By contrast, A 
should be different (see [65]). They are extracted from a s {Mz) by using b- and c-quarks thresholds following to [66]. 
The values of A obtained by this procedure and used hereafter in all the fits are: 

Alo(/ = 5) = 80.80 [MeV] , A LO (/ = 4) = 111.8 [McV] , A LO (/ = 3) = 136.8 [MeV] , 

Ams(/ = 5) = 195.7 [MeV] , K m {f = 4) = 284.0 [MeV] , A M g(/ = 3) = 347.2 [McV] , ^> 

obtained from ZEUS result a s (M z ) = 0.1166 (see [7]). 

Table II contains the results of separate fits to HI and ZEUS data with a low Q 2 cut, Q 2 ut , that increases step by 
step. We observe that the agreement between theory and experiment improves when increasing the value of Q 2 ut . F° r 
Q 2 > 2.5 GeV 2 the agreement is good (see Tables II and III). 

Note that the separated fits of HI and ZEUS data lead to purely comparable values of the parameters Q 2 ,, x , Aq , 
A t2 . Thus we may fit to the combined data set. The results of such combined fits can be found in the last rows of 
Table II and Table III. 

Looking carefully on that Tables, we arrive to the following conclusions: 

• In the leading twist approximation the preferred number of flavors / is four. 

• The value of the quark distribution does not depend on the specific Q 2 ut values within the limits of experimental 
errors. The magnitude of the gluon density and Qq decrease slowly with decreasing <2 2 ut . 

• A strong reduction of the magnitude of the gluon density is observed when NLO corrections are included. 

The suppression of the gluon density rise with Q 2 at NLO in comparison with the LO prediction is well-known 
effect [59, 67] but in addition we also observe a strong reduction of the gluon magnitude at Qq. 

At least partialy, this effect can be explained based on the GRV-like point of view [31-34], where at low Q 2 
values there are only valence quarks and all other types of partons are generated in the Q 2 -evolution. Thus, the 
slowe rise with Q 2 when NLO corrections are included directly implies a reduction of the magnitude at a given 

Ql 

It should be mentioned that a similar relative reduction of gluon normalization is obtained in the analyses 
[35, 68], when the ln(l/a;) resummation was included. Thus, the correct incorporation of NLO terms has a 
similar tendency. 

• The fitted Qq values are essentially higher at NLO: Qq <~ 0.5 -j- 0.6 GeV 2 , in comparison with LO fits, where 
Qq <~ 0.3 -j- 0.4 GeV 2 , and comparable to those obtained earlier in [44]. 

Partialy, the effect can be explained by different A values at LO and NLO approximations. Note, however, that 
the ratio A^jg/A 2 <~ 6.4 and, thus, the Q 2 dependence of F 2 data itself should be important in the definition 
ofQ 2 . 

Considering Tables II and III and Fig. 1 we find good agreement with data only at Q 2 > 2.5 GeV 2 . The situation is 
little bit worse than it was before in [44] , mainly due to the strong improvement of experimental data. To expand the 
range of applicability of our analysis to Q 2 < 2.5 GeV 2 wc add to our fits HT corrections presented in the previous 
sections. 

Let's consider both types of estimations of the HT corrections separately. 
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B. BFKL- motivated estimations for twist- four operators 

Tables IV- VI and Fig. 4 contain the results of the fits to HI and ZEUS data using Eqs. (16), (50) and (53) at LO 
and (25), (50) and (53) at NLO. 

The results demonstrate good agreement between the theoretical predictions having BFKL-like twist-four term and 
experimental data of the HI [3, 4] and ZEUS [11, 12] collaborations. 

The fits of HI [3, 4] and ZEUS [11, 12] data demonstrate a strong improvement of the agreement between theory 
and experiment (see Fig. 4), essentially at LO and in the case / = 4. 

The values of parameters in the twist-two terms do not change drastically. Qq rises 100 MeV and 150 MeV at LO 
and NLO, respectively. The gluon density in the twist-two term rises essentially and the quark distribution decreases 
slowly. The changes are compensated by a negative gluon and a positive quark twist-four magnitudes, respectively. 

We found also a tiny dependence on the real value of the parameter that supports our hypothese (see Section III) 
about the irrelevance of the exact form for the nonsingular (at n — > 1) terms in the twist-four anomalous dimensions. 

An interesting fact is that the value of the sum A^ + 4/9A^ 4 is very close to zero. Hence, HERA data do not seem to 
support a strong increase of the twist-four terms at small x, contrary to the expectation from various BFKL-motivated 
estimations [46-49]. However a small value for the the twist-four terms has also been found in a model-dependent 
analysis [69]. 

C. Renormalon model predictions for higher twist operators 

Tables IV-VII and Figs. 4 and 5 contain the results of the fits to HI and ZEUS data using Eqs. (16) and (66) 
at LO and (25) and (66) at NLO. The results demonstrate excellent agreement between theoretical predictions and 
experimental data. The x 2 decreases very strongly. 

Consider separately the fits of data for Q 2 > 1.5 GeV 2 and Q 2 > 0.5 GeV 2 , presented in Tables IV- VI (and on Fig. 
4) and Table VII (and on Fig. 5), respectively. 

Looking carefully Tables IV- VI and Fig. 4, we arrive to the following conclusions: 

• For the data usage of / = 4 is strongly preffered. 

• The values of parameters in the twist-two terms do not change essentially. 

We see, however, for HI data in Table IV and for combibed data in Table VI some rise of gluon terms when higher 
twist terms are incorporated. The rise exists for both the LO and NLO approximations and it is compensated 
by negative gluon twist-four magnitude. The twist-six gluon magnitude has different signs (it is negative and 
positive at LO and NLO approximations, respectively) but the combination of the higher twist terms gives 
negative contribution for the gluon case. 

Note that the phenomenon is similar to one observed for BFKL-motivated twist-four corrections (see previous 
subsection) and can be considered as quite general property of the HT corrections. 

• For the ZEUS data in Table V the influence of the higher twist terms is not so important. 

• In contrary to the gluon case, the higher twist corrections for the quark density are mostly positive that leads 
to different small- a; asymptotics of gluon and quark distributions at low Q 2 values, observed recently at HERA 
experiments [70] (see a detailed discussions in the subsection F). 

• The fitted value of Q 2 tends to be little higher (at LO Ql ~ 0.5 GeV 2 and at NLO Q 2 , - 0.7 0.8 GeV 2 ) when 
the twist-four corrections have been added. It is in agreement with the results when BFKL-motivated twist-four 
corrections have been considered (see the previous sunbsection) . The incorporation of twist-six terms returns 
the Qq values to the ones, obtained in the twist-two approximation. 

Looking carefully Table VII and Fig. 5, we see full support of above results: the agreement with experimental data 
improves drastically, essentially for 0.5 GeV 2 < Q2 < 2.5 GeV 2 . We should note, however, about following excepting 
features: 

• Usage of / = 3 is preffered, that is natural choice at low Q 2 values. 

• The twist-six corrections are important to stabilize the HT contributions and, thus, the results of Table VII are 
comparable with ones in Tables IV- VI only when the twist-six corrections taken into account. 
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D. Leading and higher twist approximations for the derivative dF2/d\nQ 2 

The results for the derivative dF 2 /d\nQ 2 are shown on Fig. 3 and Fig. 8 together with HI experimental data [2]. 

Fig. 3 contains only the leading twist theoretical predictions. As in the case of F 2 data we have very good agreement 
between our formulae and experimental data at Q 2 > 3 GeV 2 . 

When we added the HT corrections, the theoretical results begin to be in agreement with experiment also at 
Q 2 < 3 GeV 2 (see Fig. 8), especially when we used the results of F 2 data fits at Q 2 > 0.5 GeV 2 . The corresponding 
results for dF 2 j 'd\nQ 2 are shown as the dashed curve for the NLO and as the dash-dotted curve for the LO fits. 

Both curves are hardly distinguished from each other. It mecns, that in this kinematical region of small x the order 
of perturbation theory inside the leading twist does not matter. The importance has the number of twists taking into 
account. 

Note that the HT corrections to F 2 and dF 2 /dhiQ 2 structure functions are opposite in sign that demonstrates the 
importance, respectively, the quark density and gluon one for the functions (see also the following subsection and 
discussions therein). The fact is in full agreement with results of Section VI. 

Thus, our quite simple formulas obtained in the generalized DAS approach are very convenient also to the study 
the derivative dF 2 /d\nQ 2 , which is very important to extract gluon density and the longitudinal Fl or the ration 
R = a L /a T (see [71, 72] and [73-75], respectively). 

E. Effective slope Af£(a:, Q 2 ) 

The results for the slope X c p(x, Q 2 ) are shown on Figs. 2, 6 and 7 together with HI and ZEUS experimental data 
[1, 8, 76, 77]. 

At Figs. 2 and 6 we see very good agreement between theory and experiment as with and without consideration 
of the HT corrections. Note that the asymptotic approximation does not work so well because at large Q 2 values, i.e. 
at its range of applicability, there arc experimental data only at quite large x values: x > 10~ 3 . 

Since the logarithmic x derivative is compatible with independence of Q 2 , HI and ZEUS have both fitted their data 
on the proton structure function to the form F 2 — c(Q 2 ) x~ x (® K Fig. 7 shows recent HI and ZEUS fits [1, 8, 9, 76] 
for \{Q 2 ). Some of them arc preliminary only and extracted from Fig. 14 of the recent review [77]. 

The experimental data shows a rise of the slope \°p 2 (x, Q 2 ) from the value ~ 0.1 at Q 2 < 1 GeV 2 (so-called "soft 

pomeron range") to the value ~ 0.3 -j- 0.4 at Q 2 > 100 GeV 2 (so-called "hard pomeron range") and c ~ 0.18 is 
consistent with being constant. 

In our opinion, the strong Q 2 dependence of the slope \ e p 2 (x, Q 2 ) was observed firstly in Ref. [78], where fits of 
experimental data have been performed for the Rcgge-like PD form. At high Q 2 side, the slope value is close to 
LO BFKL prediction (Xf 2 (x,Q 2 ) ~ 0.3 + 0.4), at smaller Q 2 values \f 2 (x,Q 2 ) ~ 0.2, that is close to model with 
Pomeron interactions [79] and to NLO BFKL predictions [63] based on non-MS-likc rcnormalization schemes and 
BLM resummation of large values of NLO corrections calculated recently in [80, 81] (see also [82, 83] ). At low Q 2 
the slope value coincides with Donnachie-Landshoff model, where X c p 2 (x, Q 2 ) <~ 0.1. 

In a sence, the shape of the slope Xp 2 (x, Q 2 ) is in contrast with Regge asymptotics, where the corresponding slopes 
should be Q 2 -indcpcndcnt. Note, however, that this Q 2 -dependence can be described in phenomenological Reggc-likc 
models [84-87]. There are also attempts (see [39, 40]) to recove the slope shape in the Regge-like form of parton 
distributions considering the small x asymptotics of DGLAP equation. A quite natural explanation of the rise is given 
in the generalized DAS approximation as it was shown in [44]. 

Quite recently the HI 96/97 data [1] (black circles on Fig. 7) has been analysed in [53], where good agreement has 
been found between data and theoretical predictions based on generalized DAS approach. For example, the rise can 
be described as lnlnQ 2 , i.e. in pure perturbative QCD. Incorparation of HT corrections gives a possibility to extend 
the agreement to new preliminaty HI and ZEUS data for quite low Q 2 values (see dashed curve and the preliminary 
data near Q 2 ~ 1 GeV 2 on Fig. 7). 

F. Parton distributions 

The results for the quark and gluon densities are shown on Fig. 9 together with the NLO QCD predictions of 
A02NLO [88], represented by dots. 

As it was noted already in Sections VII and IX, there is very strong difference between the twist-two and total 
parton distributions. In the case of the twist-two parton densities f^ 2 (x,Q 2 ) the higher- twist corrections contribute 
to the Wilson coeffcient functions, i.e. (in MS-like factorization scheme used here) to the relation between the parton 
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distributions f^ 2 (x,Q 2 ) and F 2 . Then, the higher- twist terms give additional power- like corrections to the relation 
and, thus, change it. 

Contrary to this, in the case of the total parton densities f a (x,Q 2 ) the coefficient functions are pure twist-two 
ones, i.e. the relation between the parton distributions f a (x,Q 2 ) and the structure function F 2 taken in the standard 
MS- like way. Thus, in the case the higher- twist corrections are responsable for the difference between the twist-two 
parton distributions f^ 2 (x,Q 2 ) and the full ones f a (x,Q 2 ). 

For the quark density the difference between twist-two distributions and total densities are not very strong. In 
Fig. 9 one can see good agreement between quark distributions obtained in the different approximations. For the 
renormalon higher-twist corrections, our results very close to obtained by ZEUS Collaboration in [70]. 

At high Q 2 values there is also good agreement between gluon distributions obtained in the different approximations. 
For small Q 2 values, in the renormalon model our total gluon density is strictly less then twist-two one: for example, 
at Q 2 = 2 GeV 2 the ration fc(x, Q 2 )/ Jq 2 (x, Q 2 ) < 1/3. Neverseless, there is a disagreement here between our results 
and the recent one from ZEUS Collaboration (see [70]) in the range of small Q 2 values: our total gluon density is 
higher essentially the ZEUS one. A similar disagreement exist between our total gluon distribution and Alckhin one 
[88] (see Fig. 9). Thus, the deviation between our twist- two and total gluon distributions is strong but less to have 
agreement with experimental data. 

In our opinion, most part of the difference comes from neglection of valent quark part fv(x, Q 2 ) in our article. The 
neglection is a quite standart tool at small x range (and quite large Q 2 values, where parton model is applicable), 
because f v (x, Q 2 ) ~ x Xv with Ay ~ 0.3 -e- 0.5. 

At low Q 2 values, however, the ignoring the valent and nonsinglct quark distributions cannot be the correct 
approximation, because here the singlet parton distributions (at least the gluon density) start to fall when x — > 0. 
Moreover, at higher orders of perturbation theory strong double-logarithmic terms contributtc to the valent and 
nonsinglet quark distributions. The contributions can be evaluated in the framework of BFKL-like approach and they 
can lead to essential decreasing of the Xy value at low Q 2 values (see [89, 90] and discussion therein). 

Thus, in our model gluon density at small Q 2 values includes effectively a contribution of the valent quark distri- 
butions and, thus, is large essentially to compare with ZEUS and Alekhin predictions from [70] and [88], respectively. 

Note that the absence of the valent quarks can be partially responsable for some disagreement between theory and 
experiment for the derivation dF 2 /d\nQ 2 (see Figs. 3 and 8 and discussions in subsection D), which is depended 
strongly on gluon density. 

We plan to return to the study the problem and to incorporate the valent quark densities in our future investigations. 

XI. CONCLUSIONS 

In generalized DAS approximation we have incorporated HT corrections for semi-analytical solution of DGLAP 
equation obtained earlier in [44] at LO and NLO levels in the leading twist approximation for the flat initial condition. 

The HT corrections have been added in two models, the so-called BFKL-like one and the renormalon one. In 
both models the HT terms lead to improvement of the agreement with new precise experimental data of HI and 
ZEUS Collaborations. The elements of the renormalon model, however, are essentially better defined and the model 
describes experimental data much better, especially at very low Q 2 values (Q 2 > 0.5 GeV 2 ). 

After verification of all uncalculable parameters in our formulae from the fits of F 2 data we apply our approach to 
compare with HI data for the derivative dF 2 /d\nQ 2 7 with HI and ZEUS data for the effective slope X e p 2 (x, Q 2 ) data 
and with experimental predictions for the parton distributions. 

We have found rather good agreement with data for the effective slope \ c p 2 {x, Q 2 ) and for the derivative dF 2 /d\nQ 2 
and also with experimental predictions for the quark distribution, but have some disagreement with other results for 
gluon densities at low Q 2 values (see the subsection F in the previous section and discussions therein), that needs an 
additional investigations. 

As next steps we plan to add at low Q 2 values to our analysis some phenomenological models of coupling constant. 
We hope to apply the Shirkov-Solovtsov analitization [91, 92] and a "freezing" procedure (see, for example, Ref. [93] 
and discussions therein). 

Moreover we plan also to add to our initial conditions (15) corrections ~ ln(l/x) and <~ ln 2 (l/x) obeying to 
Froassart restriction by analogy with consideration of these corrections in the Regge-like small-x asymptotics of 
parton distributions done earlier in [94, 95]. 

Addition of HT terms should be important also for high-energy cosmic rays, where they can lead to quite important 
shadowing corrections for cross-sections of neutrino-proton scattering studyed in DAS approach in [96]. The subject 
will be considered in forthcoming article. 

We are considering also to extend the application of the higher twist corrections for the longitudinal srtucture 
function Fl- The consideration of Fl should be very important essentially at low Q 2 values, where Fl should go to 
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zero when Q 2 — ► [97-100] at low x values based on ^-factorization procedure [101-104]). 

In the QCD improved parton model the LO results for Fl ~ a s {Q 2 ) an d, thus, do not lead to zero values to the 
longitudinal srtucture function. Moreover, the NLO corrections to Fl are large and negative at low x values (see 
[105-112] and, thus, give large negative contributions at low Q 2 range [62, 113-115]. Thus, they can lead to the 
negative values for Fl [29, 62] of perturbation theory and needs a resummation of large corrections at low Q 2 values. 
Based on Grunberg approach [116, 117], the resummation leads to recovering well-know Callan-Gross relation Fl = 
at asymptotics x — > (see [62]). 

Thus, there are a quite conserval results for Fl at low x and Q 2 values. The incorporation of the higher-twist 
corrections, which can be very important namely in the case of the longitudinal srtucture function (see recent study 
[69] and discussions therein), should give an additional important information about Fl structure at low x and Q 2 
values. Moreover, the measurement of F L should become possible in nearest future (see discussions in Section VII of 
[15]) with the proposed updates to the HERA machine, which will yield very large integrated luminosity. Note that 
some precise preliminary results for Fl can be found already in the recent review [77] and we plan to study of them 
in nearest future. 
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XII. APPENDIX A 



The twist-four and twist-six contributions in the framework of IR-renormalon model have been calculated in [118] 
(for nonsinglct case) and in [50] (for singlet one). As we already noted in Section II, we neglect the nonsinglct 
component in our analysis. The higher-twist corrections to singlet case contain sum of nonsinglet and singlet higher- 
twist results. However, we have interest only to n — > 1 asymptotics of the corrections, where nonsinglet part of 
higher-twist corrections are neglected because exact Bjorken sum rule. 

The singlet part of higher-twist corrections may be presented in the following form 



M?(n,Q 2 ) = M a (n,Q 2 ) 



1 + — u t4 

+ Q 2 Ma 




(Al) 



The quark contributions /xl 4 (n, ln(A/Q 2 )) and ^ 6 (n, ln(A/Q 2 )) [50] may be transformed to n-space 



H T a m [ n, In 



A 



SC F T R f 

Pi 



B™{n) + bl m {n) In 



(A2) 
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where 



B T q \n) 



16 1 22 1 11 1 25 1 74 1 



15 (n- l) 2 225 n-l 3 n+1 9n + 2 75 n + 4 
12 6 4 12 

~ (n + l) 2 + (n + 2) 2 + (n + l) 3 + (n + 2) 3 ' ( ' 

h r Hn) = + _A_ _ 2 1J_ ? ^_ (A4) 

9 1 ; 15n-l n+1 'in + 2 5n + 4 (n + l) 2 (n + 2) 2 ' V j 

„ Tfi . . 128 1 572 1 52 1 32 1 16 1 

B l» = ~ 777: 77 772 - 7777^7-^ + 7777^ + 7777"^ + 



105(n-l) 2 11025n-l 75 n + 1 9n + 2 3n + 3 
724 1 452 1 16 1 16 



75 n + 4 3675n + 6 5 (n + 1) 2 (n + 3) 2 ' 



„( )^ - - 2 + n + n 2 



B ^ = - isfc (j + ii) + ° {5) • 6 ° (B) = isb + 0{5) ' 

with 5 = n — 1. 

XIII. APPENDIX B 

We present here the detailed analysis [137] of the method of replacing the convolution of two functions by a simple 
product at small x. We restrict ourselves to the accuracy O(z). Some earlier presentations can be found in [54, 55] 
(here the accuracy 0(z 2 ) has been considered) and in [44]. 

Let us to consider the set of PD with different forms: 

(I) Regge-like form f R (z) = Z - s J{z), 
(II) Logarithmicdike form /l(z) = z~ s \n(l/z)f(z), 

(III) Bcsseldike form //(z) = z~ s d\n(l/ zf^h (^2\J d\n(l / z)^ j(z) with definition (22) of the I k function, 
where f(z) and its derivative f'(z) = df(z)/dz are smooth at z — and both are equal to zero at z = 1: 

7(i) = 7(i) = . 



(A5) 



6 x__6£^ 8_L 8_L_ _8 24 1 16 1 

9 ^~105n-l 5n+l 3n + 2 + n + 3 5n + 4 + 35n + 6' ^ ' 

The gluon contributions /x^ 4 (n, ln(A/Q 2 )) and "g 6 (n, ln(A/Q 2 )) may be estimated [50] as 

uT/>,ln(A/Q 2 )) = u™(n,ln(A/Q 2 ))/ 7 ^(n) , (A7) 

where [119] 

7 £>) = -w F 



(n - l)n{n+ 1) 

is the leading contribution to the gluon-quark anomalous dimension. 

We have the interest to the asymptotics n — > 1, where the above values may be represented as 

n r4, n 4 ( 1 11 1 2291\ (n ,„ Lr4 . . 2 /l 139\ (n ,„ 

= - 15 U + 120* - 3600 j + ° W ' = 15 120 j + ° W < 

n r 6 , n 32 ( 1 143 1 870637 \ (n ,„ , r6 . . 16 /l 3217\ (n ,„ 
^<"> = + ^* ~ lini^ 6 » = - 105 U- 3360 (A8) 
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(1) At the beginning we consider the basic integral with integer integer nonnegative n values: 

J$in,z) = z n ®U{z)= [ ^-y n fi(-) , i = R,L,I. 
(a) Regge-like case. Expanding f(z) near /(0), we have 



M(n,z) = z- S 



J Z 



/(0) + -/<')(0) + 



= z 



— z 



1 



n + S 
1 



f(0) + O(z) 
1 



/(0) + 



n + 6-1 



/ (1) (0) 



1 



1 



n + S 

Using the power-like large x asymptotics 

f(z) ~ (1 — z) v when z 
the second term on the r.h.s. of Eq. (Bl) can be summed: 



A;! n + 5 — k 



/«(0) + .. 



J ifl( n > z ) = z " 



1 



n + S 



/(0) + 0(z) 



+ T{l + v-n-5) I[ ' ■ 



Consider particular cases n > 1 and n = separately: 

(al) If n > 1 then the second term in the r.h.s. of (B3) is negligible and we have 



M{n,z) = z- & 



1 



m+oiz 1 - 6 ) 



i 



f R {z) + 0{z 1 - & ) . 



(a2) If n = 0, the r.h.s. of (B3) can be rewritten as follows 



j£>(o,*) = *- 



T(-5)T{1 + V ) 
T{l + v-5) 



(Bl) 



(B2) 



(B3) 



(B4) 



/(0) 



= ^ 1 (^)/ fl (^) + 0(^- 6 ) 



where 



r(i-f)r(i + i/) 4 
r(i + i/-j) 



(B5) 



(B6) 



i.e. there is the correlation between small x and large x asymptotics of parton distributions (see [120-122]. Note that 
the value S^ 1 (z) is finite at the limit 5 — > 0: 



lim — — — — - 

5^0 5 R (z) 



= In - 



[*(1 + i/) - *(1)] = ln( - ) -p(i/) , 



(B7) 



where the Ricmannian ^-function is the logarithmic derivation of the T-function. 

Remember that the large x asymptotics are different in quark and gluon cases, the values v q w 3 and vq w 4 are 
coming from quark counting rules, what lead to p{v q ) « 11/6 and p(^g) ~ 25/12. 
(b) Logarithmic-like case. Using the simple relation 



z- 5 \n{l/z) = al(z- 5 )/d6 



we immediately obtain: 
(bl) n > 1 case: 



jS(n,z) = z- 5 ln(l/z) 
1 



1 - 



n + S 
1 

n + S 



1 



n + (5 V (n + (5)ln(l/z 

f L {z) + 0{z l - s ) 



f(0) + O(z) 



1 



(n + (5) ln(l/z) 
1 



ln(l/z) 



(B8) 
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(b2) n = case: 



jW(0,z)^z- 5 ln(l/z) 



<51n(l/z) 



f(z) + 0(z) 



= <5 Z 1 (z)/ L (z) + 0(z 1 - 5 ), 



(B9) 



where 



z A d / z 
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The value <5 i 1 (z) is also finite at the limit <5 — > 0: 



lim ; . . = — In 



^0(5 L (z) 2 V z / 21n(l/z) 

= imn L_ 

2 V z / 21n(l/z) 



^ ([*(1 + v) - *(1)] 2 - + „)- *'(1)]) 
(p(i/) 2 -j/(i/)) , 



(Bll) 



where the ^'-function is the derivation of the ^-function, p'{y q ) ~ —49/36 and p'{vg) ~ —205/144 are coming from 
quark counting rules. 

(c) Bessel-like case. Representing Bessel function in the form 



z~ 5 (B12) 



and repeating the above analysis, we have 
(cl) in the n > 1 case: 



fcSj in the n = case: 



where 



4V M 



ln(l/z) 



(B13) 



(B14) 
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dS 
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5l{z) d\n{l/z) k/2 I k ^dln(l/z) 
The value 6j (z) is also finite at the limit S — > 0: 



4 2, d — 
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(B16) 



I k \2^d\n(l/z) 

where the r.h.s. of (B16) is obtained from the expansion of the modified Bessel functions at z — > 0. 
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Note that we can represented the Eqs. (B5), (B9) and (B14) formally as follows 

1 



S- 1 f B (z) 



6b(z) 



f B (z) (B = R,L,I), 



(B17) 



which has been used in Sections III and V. 



(2) Since the HT coefficient functions B T q 4 ' 6 (n) contain the terms ~ l/(n — 1) (sec Eqs. (A3) and (A5)), we should 
consider also the second basic integral with integer nonnegative n values: 



jf}{n,z) = z n \n{l/z)®f i {z) 



1 dy 



y n HVv)h 



, i = R,L,I . 



It is easy to demonstrate that 



that symplyfies essentially the consideration of J$?(n,z). 

(a) Regge-like case. Repeating the analysis of the subsection (la), we obtain easy that 



+ z' 
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Consider particular cases n > 1 and n = separately: 

faij If n > 1 then the second term in the r.h.s. of (B3) is negligible and we have 



jf) i in,z)=z- s J ^ w m + Oiz 1 -*) 
(a2) If n = 0, the r.h.s. of (B18) can be rewritten as follows 



1 
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dd 5 R (z) 



Note that the value 5 R (z) is finite at the limit 5^0: 

1 1 



lim „ 

5 2 R {z) 



1 



lim 



(B20) 



(B21) 



(B22) 



where the value of lim^o (1/ (5jj(2)) is given in (B7). 

(b) Logarithmic-like case. Following to the subsection (lb), we obtain: 

(bl) n> 1 case: 



/z) 



(B23) 



(b2) n — case: 



4i(0,z) = &- L 2 (z) f L {z) + 0{z'- 5 ) 



(B24) 



29 



where 

1 

and the value of 1/5l{z) is given in (BIO). 

The value 5^ 1 (z) is also finite at the limit S — > 0: 



d 1 



(B25) 
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<5^0 ^(2) 



mil 



5^) a -^))-3E^ 



(pM 3 - 3^(^(1/)+ p"M) 



(B26) 



where the ^"-function is the second derivation of the ^-function, p"{v q ) ~ 251/108 and p"{vc) ~ 2035/865 are 
coming from quark counting rules. 

(c) Bessel-like case. Following to the subsection (lc), we obtain: 
(cl) in the n > 1 case: 



J$(n,z) 



{n + sy 



fi(z) + 



( c2) in the n = case: 



where 
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ln(l/z) 
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and the value of 



6 z (z) 



is given in (B15). 



The value Sj (z) is also finite at the limit 5-^0: 
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(B30) 



Note that the r.h.s. of (B16) is obtained from the expansion of the modified Bessel functions at z — > 0. 
Note that we can represented the Eqs. (B20), (B24) and (B28) formally as follows 



6 ~ 2 fs(z) = T2~r~\ /*(*) (B = R, L, I), 

°B\ Z ) 



(B31) 



which has been used in Section V. 
(3) Consider the Mellin integral 



h(z) = K{z)®f{z)= / ^ K(y) f 



f 

J z 



dy 



and define the moments of the kernel K(y) in the following form 

K n = f dy y n - 2 K{y) . 
Jo 



In analogy with part (1) we have for the Regge-like case: 
I S ,r(z) = z- 5 f dy y 5 - 1 K(y) 
= z- 5 \k 1+s f(0)+O(z) 



/( fe )(0) + ... 



N 1+s (x)f(0) + N s (z)f^(0) + . . . + 1 N 1+s - k (z)f {k) (0) + ■■■ 



(B32) 



30 



where 



N v (z) = [ dy y^- 2 K(zy) . 
Jo 



The case Ki + g = l/(n + 5) corresponds to K(y) = y n and has been already considered in part (1). In the more 
general cases (for example, Ki + g = ^(1 + 6) + 7) we can represent the "moment" Ki + g as a series of the sort 

E m =i !/(« + « + "»). 

So, for the initial integral at small x we get the simple equation: 

I StR (z) = z- s K R , 1+5 f(z)+0(z 1 - s ) = K 1+5 fnW + Oiz 1 - 6 ) , (B33) 

where the coefficient Kr^+s coincides with the one K\ + $ in the case if K n does not contain the term l/(n — 1). The 
coefficient Kr.i + $ contain the term S] i 1 (z) if the term l/(n — 1) contributed to K n . So, the function Kr^ + s can be 
represented in the form: 

Kr,i+s = K 1+5r ( z ) . (B34) 
Repeating the analysis of the subparts (b) and (c), one easily obtains 

I SiL {n,z) = K L , 1+s f L {z) + o(^—^ (B35) 



where 

K L ,i+s = Ki+5 L (z) i (B37) 
Ki,i+5 = K 1+Sl ( z ) . (B38) 

Thus, in the non-singular case (i.e. in the case when K n does not contain the term l/(n— 1)) the results of tans 
formation of the Mcllin convolution to usual products depend only on the 5 value but not on the concrete shape of 
parton distribution. The presence of the term l/(n — 1) in K n leads to the results depending on numerical value of 
5. If S is large (more precisely, if z~ b 3> const), the presence of the term l/(n — 1) in K n leads to the term 1/5 in 
the functions -ft^.i+cS (i = R,L,I) (because the term z 5 is negligible in expressions for l/Si) and the results do not 
also depend on the concrete shape of parton distribution. If 6 is small (i.e., if z~ s w 1 + <51n(l/z), that depends on 
concrete z values, of course), then the subasymptotic of parton distribution starts to play and the function -ft^i+a 
(i = R,L,I) contains the term 1/di, which is determined by the both: asymptotics and subasymptotics of parton 
distributions. 
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TABLE II: The result of the LO and NLO fits to HI (1996/97) [2] and ZEUS (1996/97) [7] data for different low Q 2 cuts. In the 
fits / is fixed to 4 flavors. 



\c£ [UC V j ^ 






vo IV 5 v J 


X / rt-v.p. 


LO(Hl 96/97 [2]) 
1.5 
2.0 
2.5 

o.O 


0.797±.022 
0.819±.022 
0.869±.024 
n non-1- nos 


0.791±.026 
0.781±.026 
0.754±.027 

U. ( oozt.Uzy 


0.304+.005 
0.309+.005 
0.319+.005 

U.OOZzt.UUD 


181/101 
139/98 
88/90 

Dl / ol 


LO(ZEUS 96/97 [7]) 
2.7 
3.5 


0.918±.031 
893+ 034 


0.754±.040 
780+ 049 


0.317+.005 
31 c i+ 006 


80/116 
76/111 


NLO(Hl 96/97 [2]) 
1.5 
2.0 
2.5 
3.5 


-.013±.015 
0.003±.015 
0.042±.017 
0.082±.020 


0.893+.028 
0.882+.028 
0.850+.029 
0.824+.032 


0.494+.009 
0.505+.009 
0.526+.010 
0. 554+. 012 


201/101 
153/98 
95/90 
63/81 


NLO(ZEUS 96/97 [7]) 
2.7 
3.5 


0.061±.023 
0.044±.025 


0.844+.044 
0.871+.046 


0.523+.011 
0.520+.012 


82/116 
78/111 


NLO(Hl[2] + ZEUS[7]) 
1.5 (r z = .963) 
2.0 (r z = .964) 
2.5 (r z = .963) 
3.5 (r z = .962) 


0.010±.013 
0.021±.013 
0.046±.013 
0.063±.015 


0.873+.024 
0.864+.024 
0.839+.024 
0.829+.026 


0.506+.007 
0.512+.007 
0.524+.008 
0.537+.008 


286/217 [204/101, 82/116] 
233/214 [154/98, 79/116] 
171/206 [95/90, 76/116] 
140/192 [66/81, 74/111] 



TABLE III: The result of the LO and NLO fits to HI [2-6] and ZEUS [7-14] data for different low Q 2 cuts and different /. 



q' 2 


[GeV^] 


> 






Aq 


A q 


Q'o [GcV 2 ] 


X*/n.o.p. 


LO (/ = 


3) 


















0.5 


(rm 




.933, rz 




.955) 


1.216+.015 


1.153+.015 


0.306+.003 


1163/667 


[488/292, 675/375] 


1.0 


(rm 




.939, r z 




.966) 


1.424+.023 


0.977+.023 


0.313+.003 


854/631 


[389/279, 465/352] 


1.5 


(rHi 




.946, rz 




.969) 


1.472+.024 


0.950+.023 


0.317+.003 


775/614 


[348/267, 427/347] 


2.0 


(rm 




.953, r z 




.971) 


1.527+.025 


0.923+.023 


0.323+.003 


673/591 


[273/252, 400/339] 


2.5 


(r H i 




.958, r z 




.971) 


1.589+.026 


0.890+.024 


0.330+.003 


580/573 


[193/236, 387/337] 


3.5 


(rHi 




.963, r z 




.971) 


1.655+.030 


0.866+.026 


0.339+.004 


501/532 


[142/210, 359/322] 


LO 


(/ = 


4) 


















0.5 


(r H i 




.934, r z 




.957) 


0.641+.010 


0.937+.012 


0.295+.003 


1090/667 


[455/292, 635/375] 


1.0 


(mi 




.940, rz 




.966) 


0.755+.015 


0.821+.019 


0.301+.003 


826/631 


[373/279, 453/352] 


1.5 


(r H i 




.947, r z 




.969) 


0.784+.016 


0.801+.019 


0.304+.003 


754/614 


[335/267, 419/347] 


2.0 


(rHi 




.953, r z 




.971) 


0.817+.017 


0.780+.019 


0.310+.003 


659/591 


[264/252, 395/339] 


2.5 


(rm 




.958, r z 




.971) 


0.855+.017 


0. 754+. 020 


0.316+.003 


570/573 


[188/236, 382/337] 


3.5 


{rm 




.963, r z 




.971) 


0.892+.020 


0.737+.021 


0.325+.004 


495/532 


[140/210, 355/322] 


NLO (/ 


= 3) 
















0.5 


{rm 




.929, r z 




.951) 


-.094+.009 


1.358+.015 


0.515+.006 


1406/667 


[599/292, 807/375] 


1.0 


{rm 




.936, r z 




.965) 


0.072+.014 


1.114+.024 


0.526+.006 


966/631 


[455/279, 511/352] 


1.5 


{rm 




.944, r z 




.968) 


0.109+.015 


1.078+.025 


0.535+.006 


863/614 


[403/267, 460/347] 


2.0 


{rm 




.952, r z 




.971) 


0.151+.016 


1.045+.025 


0.548+.006 


735/591 


[311/252, 424/339] 


2.5 


{rm 




.958, r z 




.970) 


0.198+.016 


1.006+.025 


0.564+.006 


620/573 


[213/236, 407/337] 


3.5 


{rm 




.963, r z 




.971) 


0.254+.019 


0.972+.027 


0.587+.007 


523/532 


[151/210, 372/322] 


NLO (/ 


= 4) 
















0.5 


{rm 




.932, rz 




.955) 


-.142+. 006 


1.087+.012 


0.478+.006 


1229/667 


[514/292, 715/375] 


1.0 


{rm 




.938, r z 




.966) 


-.042+. 011 


0.929+.021 


0.487+.006 


884/631 


[407/279, 477/352] 


1.5 


{rm 




.946, r z 




.969) 


-.020+. 011 


0.903+.021 


0.495+.006 


798/614 


[363/267, 435/347] 


2.0 


{rm 




.953, r z 




.971) 


0.006+.012 


0.877+.021 


0.506+.006 


688/591 


[282/252, 406/339] 


2.5 


{rm 




.958, r z 




.971) 


0.035+.012 


0.847+.022 


0.520+.006 


589/573 


[197/236, 392/337] 


3.5 


{rm 




.963, rz 




.972) 


0.065+.014 


0.826+.023 


0.539+.007 


505/532 


[143/210, 362/322] 
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TABLE IV: The result of the LO and NLO fits to HI (1996/97) [2] data. Power corrections included for different values of the 
parameter 6 and in the infrared renormalon case. 



HI 96/97 [2] 


AT-2 

A G 




Ag (aS 4 ) 
(a& 6 ) 


A? (aD 
(aj 6 ) 


Qo [GeV z ] 


X'/n.o.p. 


L( J (J = 4 ) 
no hr 
b = 
b = 1 
b = a 2 /2 


0.797±.022 
1.214±.060 
1.263±.070 
1.321±.072 


0.791±.026 
0.426±.054 
0.436±.051 
0.446±.049 


-.496±.062 
-.523±.065 


0.969±.127 
1.127±.142 
1.205±.148 


0.304±.005 
0.360±.009 
0.388±.022 
0.417±.023 


181/101 
124/101 
119/101 
106/101 


renorm. 


1 1 ^4- Dfif) 
J — Luuzn.uuu 

1.037±.121 


0.668±.073 


— n+ 1 71 
(0.000 fix) 
-.011±.259 
-.486±.841 


fl 9^04- D7S 
u.^ouzn.u / o 

(0.000 fix) 

-.007±.122 

0.084±.325 


0.356±.035 


54/101 


no hT 
6 = 
b = 1 
b = a 2 /2 


-.013±.015 
-.024±.017 
0.316±.047 
0.336±.045 


0.893±.028 
0.882±.029 
0.474±.056 
0.492±.053 


-.542±.065 
-.603±.067 


-.OOli.OOO 
1.219±.147 
1.362±.152 


0.494±.009 
0.473±.017 
0.600±.030 
0.635±.030 


201/101 
199/101 
133/101 
127/101 


renorm. 


0.144±.078 
0.102±.086 


0.764±.056 
0.800±.066 


-.692±.275 
(0.000 fix) 
-1.327±1.218 
0.412±.834 


0.155±.021 
(0.000 fix) 
0.310±.281 
0.063±.144 


0.576±.060 
0.548±.067 


55/101 
54/101 



TABLE V: The result of the LO and NLO fits to ZEUS (1996/97) [7] data. Power corrections included for different values of the 
parameter b and in the infrared renormalon case. 



ZEUS 96/97 [7] 


A G 


AT-2 

A i 


A h* (a& 4 ) 
(a& 6 ) 


K {<) 
K 6 ) 


Qo [GeV 2 ] 


X z /n.o.p. 


LO (/ = 4) 
no hr 
6 = 
6=1 
6 = a 2 12 


0.918±.031 
0.891±.067 
0.910±.074 
0.920±.069 


0.754±.040 
0.780±.070 
0.780±.068 
0.786±.066 


0.046±.101 
0.083±.117 


-.093±.203 
-.101±.229 
-.179±.263 


0.317±.005 
0.314±.009 
0.324±.023 
0.330±.019 


80/116 
80/116 
79/116 
78/116 


Rt4 
renorm. 


0.980±.063 
0.859±.087 


0.739±.050 
0.757±.074 


0.344±.329 
(0.000 fix) 
-2.43941.207 
-10.66±3.60 


-.1374.135 
(0.000 fix) 
1.014±.559 
4.99±1.78 


0.343±.021 
0.281±.024 


78/116 
68/116 


NLO (/ = 4) 
no hr 
6 = 
6=1 
6 = a 2 j2 


0.061±.023 
0.067±.030 
0.062±.015 
0.071±.055 


0.844±.044 
0.849±.046 
0.859±.026 
0.866±.073 


0.020±.002 
0.046±.122 


-.001±.002 
-.044±.005 
-.101±.275 


0.523±.011 
0.533±.034 
0.534±.017 
0.549±.037 


82/116 
81/116 
81/116 
80/116 


Rt4 
renorm. 


0.083±.081 
-.3294.068 


0.823±.078 
1.242±.094 


-.046±.313 
(0.000 fix) 
-1.599±.643 
-16.00842.451 


0.016±.041 
(0.000 fix) 
-.1774.173 
2.253±0.492 


0.533±.054 
0.312±.027 


81/116 
64/116 
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TABLE VI: The result of the LO and NLO fits to HI [2-6] and ZEUS [7-14] data at Q 2 > 1.5 GeV 2 . Power corrections included 
for different values of the parameter 6 and in the infrared renormalon case. 



HI [2-6] + ZEUS [7-14] 
Q 2 > 1.5 GeV 2 




A q 


AT (&T) 
(a& 6 ) 


AT (aj 4 ) 
«') 


Q'o [GeV 2 ] 


X 2 /n.o.p. 


LO (/ = 3) 

no hr (rm = -946, r z = .969) 
6 = (r H i = .953, r z = .970) 
6 = 1 {r m = -954, r z = .971) 
b = a 2 12 (rm = -954, r z = .972) 


1.472±.024 
2.083±.056 
2.164±.068 
2.224±.067 


0.950±.023 
0.513±.042 
0.528±.040 
0.546±.039 


-.623±.049 
-.617±.051 


1.275±.103 
1.422±.113 
1.431±.116 


0.317±.003 
0.372±.006 
0.401±.014 
0.421±.013 


775/614 [348/267, 427/347] 
628/614 [246/267, 382/347] 
616/614 [240/267, 376/347] 
591/614 [224/267, 367/347] 


JXT4: 

(rm = .959, r z = .973) 
renorm. 
(rm = .959, r z = .972) 


o ai ox_ nfio 
z.Uiz±.Uoz 

1.826±.100 


U.D8 ( ±.Uzy 
0.784±.050 


— .z <y±.ioo 
(0.000 fix) 
-.064±.185 
-1.245±.718 


U.ozo±.Uoo 
(0.000 fix) 

0.006±.092 
.525±.383 


A QAAX- A1 O 

u.oyu±.uiz 
0.360±.019 


rno /£? 1 A 

OUo/014 

[151/267, 352/347] 
498/614 

[149/267, 349/347] 


LO (/ = 4) 

no hr (rm = .947, r z = .969) 
6 = (mi = .953, r z = .970) 
6 = 1 (rm = .954, r z = .971) 
6 = a 2 /2 (rm = .954, r z = .972) 


0.784±.016 
1.157±.036 
1.202±.042 
1.232±.041 


0.801±.019 
0.461±.035 
0.477±.033 
0.494±.032 


-,478±.040 
-.481±.042 


0.950±.082 
1.088±.091 
1.107±.096 


0.304±.003 
0.353±.005 
0.383±.014 
0.402±.013 


754/614 [335/267, 419/347] 
625/614 [244/267, 381/347] 
612/614 [237/267, 375/347] 
586/614 [220/267, 366/347] 


(r H i = .959, r z = .973) 
renorm. 
(rm = .959, r z = .972) 


1 1 OCX flQ7 

0.990±.060 


U.OOZ±.Uz4 

0.679±.043 


— .1 / Zzh.lUU 
(0.000 fix) 
-.009±.161 
-.980±.497 


U.loz±.U4 ( 
(0.000 fix) 
-.019±.092 
0.276±.196 


A A1 O 

U.O I o±.Ulz 

0.345±.017 


trnc /£? 1 A 
OUO/014 

[153/267, 352/347] 
497/614 
[149/267, 348/347] 


"NTT /~ \ ( f* o\ 

NLO (/ = 3) 

no hr (rm = .944, r z = .968) 
6 = (rm = .945, r z = .968) 
6 = 1 (rm = .953, r z = .970) 
6 = a 2 /2 (r H i = .953, r z = .970) 


0.109±.015 
0.101±.018 
0.600±.043 
0.630±.042 


1.078±.025 
1.073±.025 
0.563±.043 
0.585±.041 


-.735±.052 
-.783±.053 


-.OOli.001 
1.655±.117 
1.771±.120 


0.535±.006 
0.527±.011 
0.661±.020 
0.691±.019 


863/614 [403/267, 460/347] 
862/614 [401/267, 461/347] 
669/614 [273/267, 396/347] 
654/614 [265/267, 389/347] 


HT4: 

(rm = .960, r z = .973) 
renorm. 
(r H i = .959, r z = .972) 


a a i ax_ nfi'? 
U.41U±.U0 1 

0.188±.080 


U.004±.U41 

0.985±.045 


(?(?nx ice 

— .ddu±.1oo 
(0.000 fix) 
-3.081±1.019 
.036±.422 


A Ofi C J_ AO /I 

U.zOD±.Uz4 
(0.000 fix) 
.957±.315 
.437±.164 


U.D4o±.Uo4 

0.530±.041 


OUO/014 

[151/267, 354/347] 
498/614 
[149/267, 349/347] 


ntt r» /" f /n 

no (rm = -946, r z = .969) 
6 = (r H i = .946, r z = .969) 
6 = 1 (r H i = .953, r z = .970) 
6 = a 2 /2 (r H i = .954, r z = .971) 


-.020±.011 
-.024±.013 
0.288±.029 
0.301±.028 


0.903±.021 
0.899±.022 
0.515±.037 
0.535±.035 


-.535±.042 
-.580±.044 


0.00±.0004 
1.205±.095 
1.311±.100 


0.495±.006 
0.488±.011 
0.602±.019 
0.631±.019 


798/614 [363/267, 435/347] 
798/614 [362/267, 436/347] 
645/614 [256/267, 389/347] 
629/614 [248/267, 381/347] 


Rt4 

(rm = .960, r z = .973) 
renorm. 
(rm = .959, r z = .973) 


0.156±.041 
0.041±.045 


0.734±.035 
0.824±.034 


-.522±.153 
(0.000 fix) 

-2.765±.968 
.939±.718 


0.141±.014 
(0.000 fix) 
.676±.240 
.252±.099 


0.579±.031 
0.493±.037 


506/614 

[151/267, 355/347] 
500/614 

[151/267, 349/347] 



37 



TABLE VII: The result of the LO and NLO fits to HI [2-6] and ZEUS [7-14] at Q 2 > 0.5 GeV 2 . Power corrections included in the 
infrared renormalon case. 



HI [2-6] + ZEUS [7-14] 
Q 2 > 0.5 GeV 2 


A G 


ATI 

A q 


a G 


a 9 


Ql [GeV 2 ] 


X 2 /n.o.p. 


LO Rt4 (/ = 3) 

(rm = .953, r z = .975) 
LO (/ = 3) 

(Ytm = 955 rv — 9741 


2.212±.050 
2.161±.055 


0.602±.027 
0.633±.029 


0.238±.019 
(0.000 fix) 
-.002±.024 
— .017±.017 


-.014±.005 
(0.000 fix) 
0.165±.024 
0.053±.010 


0.428±.008 
0.421±.010 


569/667 

[192/292, 377/375] 
553/667 

[181/292, 372/375] 


LO Rt4 (f = 4) 

(r m = .953, r z = .975) 
LO (/ = 4) 

(rm = -955, r z = .974) 


1.234±.031 
1.211±.033 


0.518±.023 
0.539±.023 


0.201±.016 
(0.000 fix) 
-.002±.020 
0.001±.010 


-.011±.003 
(0.000 fix) 
0.102±.015 
0.031±.005 


0.407±.008 
0.404±.009 


573/667 

[193/292, 380/375] 
555/667 

[182/292, 373/375] 


NLO RtA (f = 3) 

(rm = -951, r z = .975) 
NLO (/ = 3) 

(rm = .956, r z = .975) 


1.014±.057 
0.617±.058 


0.521±.034 
0.742±.038 


0.632±.048 
(0.000 fix) 
-.129±.102 
-.203±.053 


0.188±.026 
(0.000 fix) 
0.224±.022 
0.061±.010 


0.956±.031 
0.746±.030 


621/667 

[220/292, 401/375] 
562/667 

[182/292, 380/375] 


NLO Rt4 (f = 4) 

(rm = .950, r z = .975) 
NLO (/ = 4) 

(rm = .955, r z = .974) 


0.485±.033 
0.279±.038 


0.476±.029 
0.640±.034 


0.556±.042 
(0.000 fix) 
-.143±.100 
-.044±.050 


0.071±.008 
(0.000 fix) 
0.140±.015 
0.043±.007 


0.826±.026 
0.672±.029 


617/667 

[222/292, 395/375] 
565/667 

[184/292, 381/375] 
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FIG. 1: FJ 2 {x,Q 2 ) as a function of x for different Q 2 bins. The experimental points are from HI [2-6] (open points) and 
ZEUS[7-14] (solid points). The solid, black line represents the NLO fit with x 2 / n - d - f - = 798/611 = 1.31 [A T G 2 = -.020, A T q 2 = 
.903, Ql = .495 GeV 2 ]. The long dashed, red line represents the LO fit with xV n - d - f - = 754/611 = 1.23 [Ag = .784, A T q 2 = 
.801, Ql = .304 GeV 2 ]. 
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FIG. 2: The derivative function (effective slope) A^ t2 (x,Q 2 ) = dlnFp(x, Q 2 )/d\n (1/x) as a function of x for different 
Q 2 bins. The experimental points are from HI [2]. The outer error bars include statistical and systematical errors added in 
quadrature, while the inner error bars correspond to statistical errors only. The solid, black line represents the NLO fit with 
X 2 / n - d - f - = 798/611 = 1.31 [A r G 2 = -.020, A T q 2 = .903, Q 2 . = .495 GeV 2 ], while the long dashed, red line (hardly distinguished 
from the solid one) is the LO fit with x 2 / n - d - f - = 754/611 = 1.23 [A}? = .784, A T q 2 = .801, Ql = .304 GeV 2 ]. The dotted, blue 
line corresponds to the asymptotic expression \ a p ,T2 a (x,Q 2 ) in Eq. (43c). 
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FIG. 3: The derivative function 9i ? J 2 (a;, Q 2 )/d In Q 2 taken at fixed Q 2 and plotted as a function of x. The experimental 
points are from HI [2]. The outer error bars represent the quadratic sum of statistical and systematical errors. The inner 
error bars show the statistical error only. The solid, black line represents the NLO fit with x 2 /n.d.f. = 798/611 = 1.31 
[Ag = -.020, A T q 2 = .903, Qo = .495 GeV 2 ], while the long dashed, red line is the LO fit with x 2 / n - d - f - = 754/611 = 1.23 
[Aq = -784, A T q 2 = .801, Ql = .304 GeV 2 ]. 
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FIG. 4: F 2 (x,Q 2 ) as a function of x for different Q 2 bins. The experimental points are the same as on Figure 1. The solid, 
black line represents the NLO fit alone with x 2 / n - d - f - = 798/611 = 1.31 [A T J = -.020, A T q 2 = .903, Q\ = .495 GeV 2 ]. The 
dash-dotted, red curve represents the BFKL-motivated estimations for higher twist contribution to F2(x,Q 2 ) with the value 
of the parameter b = a 2 /2. The corresponding xV n - d - f - = 629/609 = 1.03 [A}? = .301, A\ 2 = .535, Q% = .631 GeV 2 and 
A T c = -.580 GeV 2 , A T q 4 = 1.311 GeV 2 ]. The dashed, blue curve is obtained from the fits at the NLO, when the renormalon 
contributions of higher-twist terms have been incorporated. The corresponding x 2 /n.d.f. = 500/607 = 0.82 [Aq — .041, A T q 2 = 
.824, Q 2 = .493 GeV 2 and ag 1 = -2.765 GeV 2 , aj 4 = .676 GeV 2 , ag 6 = .939 GeV 4 , a T q 6 = .252 GeV 4 ]. 




FIG. 5: F 2 (x,Q 2 ) as a function of x for different Q 2 bins. The experimental points are from HI [2-6] (open points) and 
ZEUS[7-14] (solid points). The solid, black line represents the NLO fit alone with xVn.d.f. = 798/611 = 1.31 L4g 2 = 
-.020, A T q 2 = .903, Ql = .495 GeV 2 ]. The dashed, blue curve is obtained from the fit at the NLO, when the renormalon 
contributions of higher-twist terms have been incorporated. The corresponding x 2 /n.d.f. = 565/660 = 0.86 L4g 2 — .279, A r2 = 
.640, Ql = .672 GeV 2 and ag 4 = -.143 GeV 2 , aj 4 = .140 GeV 2 , ag 6 = -.044 GeV 4 , a T q 6 = .043 GeV 4 ]. The dash-dotted, 
red curve (hardly distinguished from the dashed one) represents the fit at the LO together with the renormalon contributions 
of higher-twist terms. The corresponding x 2 / n -d.f. = 555/660 = 0.84 [Aq = 1.211, A T q 2 = .539, Q 2 , = .404 GeV 2 and 
ag* = -.002 GeV 2 , aj 4 = .102 GeV 2 , ag 3 = .001 GeV 4 , a T q 6 = .031 GeV 4 ]. 
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FIG. 6: The derivative function (effective slope) \f 2 = din F 2 {x,Q 2 )/dhi(l/x) as a function of x for different Q 2 bins. The 
experimental points and the solid, black line (NLO fit with x 2 / n -d.f. = 798/611 = 1.31) are the same as on Figure 2. All 
other curves are obtained from the fits, when the renormalon contributions of higher-twist terms have been incorporated. The 
dashed, red one is the same as on the Figure 5 with the corresponding x 2 / n -d.f. = 565/660 = 0.86, while the dash-dotted, blue 
line is the one from the Figure 4 with x 2 / n -d.f. = 500/607 = 0.82. The dotted, red line corresponds to the asymptotic LO 
expression Xf 2iS3 (x,Q 2 ) in Eq. (79), plotted at x 2 /n.d.f. = 573/667 = 0.87 [Aq = 1.234, A T q 2 = .518, Ql = .407 GeV 2 and 
aS 4 = .201 GeV 2 , a T q 4 = -.011 GeV 2 ]. 
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FIG. 7: The derivative function (effective slope) Xf 2 = d In F 2 (x, Q 2 )/d\n (l/x) as a function of Q 2 . The experimental points 

are those HI and ZEUS have fitted their x < 0.01 data to the form F 2 = c(Q 2 )x~ mq2) : black points - HI F 2 data [1]; blue 
squares - HI data [76] combined with NMC data [123]; blue triangles - HI data [76] combined with low Q 2 ZEUS BPT data 
[8]; open red diamonds - preliminary ZEUS slope fit 2001 [77]. The inner error bars illustarte the statistical uncertainties, 
the full error bars represent the statistical and systematic uncertanties added in quadrature. The data are compared with 
a parametrization [1] in which X(Q 2 ) = aln[Q 2 / A 2 ] grows logarithmically with Q 2 [a — .0481, A = 292 MeV], using data 

short-dashed one) 



for Q 2 > 3.5 GeV 2 (black solid straight line goes to approximated short-dashed one). The solid, black line (NLO fit with 
X 2 / n - d - f - = 798/611 = 1.31), the long-dashed, red one (NLO&Rht fit with xV n - d - f - = 565/660 = 0.86) and the dash-dotted, 
blue one (NLO&Rht fit with x 2 / n -d-f- = 500/607 = 0.82) are the same as on the previous Figure 6. The value of x was fixed 
to 10~ 3 for all curves. 




1 o" 5 1 o" 4 1 0" 3 1 o" 2 1 0" 5 1 o" 4 1 0" 3 1 o" 2 1 0" 5 1 o" 4 1 0" 3 1 o" 2 1 0" 5 1 o" 4 1 0" 3 1 o" 2 1 



X 



FIG. 8: The derivative function 8F2(x,Q 2 )/dlnQ 2 taken at fixed Q 2 and plotted as a function of a;. The experimental points 
and the solid, black line (NLO fit with x 2 / n -dT = 798/611 = 1.31) are the same as on Figure 3. The dashed and dotted, blue 
curves are obtained from the fit at the NLO, when the renormalon contributions of higher- twist terms have been incorporated. 
The dashed one is the same as on the Figure 5 with the corresponding x 2 /n.&.f. = 565/660 = 0.86, while the dotted line is 
the one from the Figure 4 with x 2 / n -d.f. = 500/607 = 0.82. The dash-dotted, red curve (hardly distinguished from the dashed 
one) is the same as as on the Figure 5 and represents the fit of data on structure function F2(x, Q 2 ) at the LO, the renormalon 
contributions of higher-twist terms included. The corresponding y^ 2 /n.A.i. — 555/660 = 0.84. 
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FIG. 9: The parton distributions f a (x,Q 2 ) as a function of x for Q 2 = 2, 10 and 100 GeV 2 compared to the NLO QCD 
predictions of A02NLO [88], represented by black dots. The solid lines represent the NLO fit alone with x 2 / n -d.f. = 798/611 = 
1.31 [Aq = -.020, Ag 2 = .903, Q% = .495 GeV 2 ]. The dash-dotted curves represent the BFKL-motivated estimation for 
the higher twist contribution with the value of the parameter b = a 2 /2. The corresponding % 2 /n.d.f. = 629/609 = 1.03 
[Aq = .301, A T q 2 = 0.535, Ql = .631 GeV 2 and A T G 4 = -.580 GeV 2 , A r q 4 = 1.311 GeV 2 ]. The dashed curves are obtained 
from the fits at the NLO, when the renormalon contributions of higher- twist terms have been incorporated. The corresponding 
X 2 /n.d.f. = 565/660 = 0.86 [A T G 2 = .279, A\ 2 = .640, Ql = .672 GeV 2 and ag 1 = -.143 GeV 2 , aj 4 = .140 GeV 2 , a£? = 
-.044 GeV 4 , aj 6 = .043 GeV 4 ]. 



